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'^ . Abstract 

Ch ■ 

Q ■ Disordered noninteracting quasiparticles that are governed by a 

CJ I Majorana-type Hamiltonian - prominent examples are dirty super- 

conductors with broken time-reversal and spin-rotation symmetry, or 
^^ ■ the fermionic representation of the 2d Ising model with fluctuating 

;<— ^ . bond strengths - are called class D. In two dimensions, weakly dis- 

OO I ordered systems of this kind may possess a metallic phase beyond the 

^^ ■ insulating phases expected for strong disorder. We show that the 2d 

metal phase emanates from the free Majorana fermion point, in the di- 
GN I rection of the RG trajectory of a perturbed WZW model. To establish 

this result, we develop a supersymmetric extension of the method of 
nonabelian bosonization. On the metallic side of the metal-insulator 
transition, the density of states becomes nonvanishing at zero energy, 
by a mechanism akin to dynamical mass generation. This feature is 
explored in a model of A^ species of disordered Dirac fermions, via the 



O . mapping on a nonlinear sigma model, which encapsulates a Z2 spin 

^ ' degree of freedom. We compute the density of states in a finite system, 

and obtain agreement with the random-matrix prediction for class D, 
K^ . in the ergodic limit. Vortex disorder, which is a relevant perturbation 

H I at the free- fermion point, changes the density of states at low energy 

and suppresses the local Z2 degree of freedom, thereby leading to a 

different symmetry class, BD. 
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1 Introduction 

As is well-known, the two-dimensional Ising model has a magnetic phase 
transition, the critical behavior of which is governed by the relativistic field 
theory of a massless Majorana fermion or, on squaring the partition function, 
a massless Dirac fermion. A theme of some debate over the last fifteen years 
has been the effect of disorder on this phase transition. Disorder, when 
introduced in the form of spatial inhomogeneities in the bond strengths, is 
known to add a random mass term to the Majorana Lagrangian. Perturbative 
renormalization group calculations by Dotsenko and Dotsenko |jl]] showed 
that randomness in the mass is a marginally irrelevant perturbation leading 
to logarithmic corrections to the pure Ising critical behavior. In particular, 
the singularity in the specific heat of the Ising model persists (in a weakened 
form) in the presence of disorder. This conclusion was confirmed by Shankar 
and Ludwig , who also calculated the effect of disorder on the moments 
of the spin-spin correlation function. 

The picture emerging from this work is that the 2d Ising model with 
weakly disordered bond strengths undergoes a phase transition controlled by 
the pure Ising fixed point, and the only effect the disorder has are logarithms. 
Although this picture came to be widely accepted, a constant challenge has 
been the work of Ziegler [^ |^, 0, who claims that the thermodynamic sin- 
gularities of the Ising model are rounded off by the disorder. Specifically, he 
argues that the vanishing density of states (DoS) of the relativistic fermion 
at zero energy becomes finite as a result of nonperturbative effects. To mo- 
tivate this scenario, he refers to a model of disordered Dirac fermions on the 
lattice, for which the DoS at -E = can be proved to be strictly positive. 

Ziegler's field-theoretic calculation imitates the standard treatment of dis- 
ordered metals, and is based on a Hubbard-Stratonovich transformation in- 
troducing a composite field Q, followed by a saddle-point approximation. 
Such a strategy is in principle not unreasonable. Indeed, in the case of disor- 
dered metals the effective theory for Q ~ ipip, where ip is the basic electron 
field, not only allows the systematic study of quantum interference correc- 
tions to metallic behavior, but also captures the nonperturbative physics of 
localization in low dimension. By analogy, one expects that in the present 
case, too, an advantage might be gained by transforming from Majorana 
fields ifj to composite fields Q ~ ipip. 

An important lesson learned from disordered metals is that, when using 
the field-theoretic formulation, we must exercise particular care to correctly 



implement the symmetries of the disordered Hamiltonian. For example, any 
inaccuracy in the treatment of time-reversal invar iance (when present) will 
spoil the weak localization effects due to the cooperon mode. What, then, 
are the symmetries of the disordered Majorana theory? The situation is 
most easily explained if we switch from two real (Majorana) fermions to the 
equivalent representation by one complex (Dirac) fermion, and is as follows. 
The first-quantized Hamiltonian for a two-dimensional Dirac fermion with 
random mass m{x), 

H=( z^^'K ~'^vfO , (1) 

\^ —lOi + 02 —m[x) J ^ ' 

has a symmetry of the "particle-hole" type: 

H = -a,H^a^, a, = (j j") . (2) 

According to the general classification of disordered single-particle systems 
d^, §, this symmetry places the Hamiltonian (|l]) in class D.Q From the 
analysis of Ref . [Q , one then infers that the supersymmetric field-theory rep- 
resentation for random-mass Dirac fermions has a global invariance under 
the orthosymplectic Lie supergroup OSp(2?T,|2r2),0if the energy vanishes and 
the disorder average of a product of n Green functions is to be calculated. 
Actually, the existence of this orthosymplectic symmetry was first pointed 
out by D. Bernard [^ in his Cargese lectures on the application of confor- 
mal field-theory techniques to two-dimensional disordered systems. (Random 



Dirac fermions have also been treated by CFT techniques in |jTO[-) 

It is clear that the existence of such a symmetry will have an important 
bearing on the low-energy physics. For one thing, it was shown in |^ that 
the Gaussian random-matrix ensemble for class D has a positive density of 
states at zero energy. (This can be understood from the fact [^ that the 
eigenvalues of a Gaussian random matrix in class D behave as a noninter- 
acting gas of harmonically confined fermions on the half-line, with Neumann 

^The equation (0) fixes an orthogonal Lie algebra in even dimension, so(2A^), which is 
denoted by D^ in Cartan's table. Hence the name "I?" . 

^Throughout this paper, supergroups such as GL(n|n) and OSp(2n|2n) are understood 
to be defined over the complex number field, C, unless specified otherwise. In places 
where this fact is of particular importance, we will switch to the notation Gljc(n\n) and 
OSpc(2n|2n). 



boundary conditions at the origin.) For another, the density of states enjoys 
the status of an "order parameter" in the field-theoretic formahsm. It is well 
known that, when the condensation of an order parameter is associated with 
a spontaneous breaking of global symmetries, there appear massless modes 
due to Goldstone's theorem. In the present context, the saddle-point value of 
the Q-field (the order parameter) breaks OSp(2n|2n) symmetry, and Gold- 
stone's theorem hence forces the existence of a supersymmetric OSp(2n|2n) 
multiplet of Goldstone modes. (As usual, these modes have a physical in- 
terpretation as diffusion-like modes. They are characteristic of class D, and 
were called the "spin-singlet D-type diffuson" in |^.) It is a perplexing fact 
that these Goldstone modes appear nowhere in Ziegler's work. Nonetheless, 
they exist, and because they do, the fate of the system in the thermody- 
namic limit cannot be decided on the basis of a plain saddle-point analysis, 
but is determined by the notoriously subtle problem of interacting Goldstone 
modes in two dimensions. Solving this problem requires the use of the renor- 
malization group. Thus we are led to ask: once we have augmented Ziegler's 
approach by incorporating the orthosymplectic symmetry of the disordered 
Majorana theory, what is the prediction for the local density of states at 
E = Ql Does it vanish in the thermodynamic limit, or does it not? 

Further motivation for the present paper comes from several issues beyond 
the density-of-states controversy. According to the general scheme of |]^, an- 
other realization of symmetry class D is by the low-energy quasiparticles 
of disordered superconductors with broken time-reversal and spin-rotation 
invariance. Physically, such a situation may occur in spin-triplet supercon- 
ductors, or in spin-singlet superconductors with spin-orbit scattering, when 
time-reversal symmetry is broken spontaneously or by a magnetic field. Al- 
ternatively, both time-reversal and spin-rotation invariance can be broken by 
randomly adding classical Heisenberg impurity spins. With neither the elec- 
tric charge nor the spin of a single quasiparticle being conserved in class D, 
the only constant of the motion is the energy. Hence, quasiparticle diffusion 
and localization in a superconductor of class D has to be probed via thermal 
transport (or transport of energy). 

The qualitative physics of class-D quasiparticles was the subject of a 
recent paper by Senthil and Fisher |]ri|. These authors drew a schematic 
phase diagram for the two-dimensional case in particular. There are three 
phases, namely thermal insulator, thermal quantum Hall fluid, and thermal 
metal, with three distinct phase boundaries, which meet in a multicritical 
point, M* . The structure of the phase diagram is determined on very basic 



grounds, by localization of all states for strong disorder (giving the insulat- 
ing phases), the renormalization group flow of a weakly coupled nonlinear 
sigma model (giving the metal), and the topological distinction between the 
insulator and the quantum Hall fluid with edge currents. 

In a more speculative attempt, Senthil and Fisher went on to try and 
match the phases of disordered 2d quasiparticles in class D with the phase di- 
agram of the 2d random-bond Ising model. Such an identification is prompted 
by the representation of the 2d Ising model by a Majorana fermion. However, 
there appeared to be a difficulty: the three stable phases of symmetry class 
D in two dimensions seem to have only two counterparts, namely the ferro- 
magnet and the paramagnet, in the Ising model. Where is the missing third 
phase? To resolve this puzzle, Senthil and Fisher suggested two alternative 
scenarios. In brief, the first one identifies the multicritical point M* with a 
point on the Nishimori line [0, while the second one interprets the missing 
phase as a spin glass phase. (The latter, however, is thought to exist only 
at zero temperature in d = 2.) Neither scenario looks convincing, which will 
motivate us to take a fresh look at the puzzle. 



The latest contribution to the subject was made by Read and Green |13 
who expounded the fact that another physical realization of the massless 
Majorana theory exists (in mean-field approximation) at the transition be- 
tween paired quantum Hall states, and in chiral p-wave superconductors at 
the transition between the topologically distinct phases of weak and strong 
pairing. Concerning the role of disorder, they suggested to extend the defi- 
nition of class D so as to include vortices. Moreover, they argued that ran- 
domly placed vortices are a strongly relevant perturbation at the free-fermion 
point. The question, however, exactly what the renormalization group flows 
to, was not answered conclusively. Based on symmetry grounds only. Read 
and Green wrote down a nonlinear sigma model which differs from the one 
we obtain for class D in that the target space in the fermion-fermion (FF) 
sector, 0(2n)/U(n), is replaced by its connected component S0(2n)/U(n). 
Thus we may ask: what happened to the Z2 degree of freedom that distin- 
guishes between 0(2n) and S0(2n)? Giving an answer to this subtle and 
yet pertinent question provides the final motivation for writing the present 
paper. 

In view of the length of the paper, we now summarize our main results. 



1.1 Overview 

For generality and better perspective, we study a family of models with 
A^ > 1 species of Dirac fermions, replacing the random mass m{x) for a single 
species by a mass matrix Mki{x) for A^ species. A perturbative renormaliza- 
tion group calculation shows that in addition to the random-mass coupling, 
a few other couplings are generated by the RG flow for N > 1. By extending 
the initial formulation of the model, we incorporate the most important one 
of these. We then use standard technology to approximate the disordered 
Dirac theory by an effective action for nonlinear fields Q : M^ ^ X. As 
expected from |^ , the target space X is a Riemannian symmetric superspace 
of type CI|DIII. Its bosonic base (or "body") is a product Mb x Mp where 
Mb = Sp(2n,M)/U(n) is a noncompact symmetric space of type CI, and 
Mp = 0{2n)/U{n) is a compact symmetric space of type Dili. Note that 
by 0{2n) we do mean the full orthogonal group, which consists of two con- 
nected components (the orthogonal matrices with determinant +1 or —1). 
Thus the present target space is distinguished by the striking feature of hav- 
ing two disjoint components, and there exists the possibility, not previously 
encountered in this context, of forming Z2 domain walls in the Q- field theory. 
The effective action for the nonlinear field Q is the logarithm of a su- 
perdeterminant, resulting from integration over the Dirac fields. One now 
wants to expand the effective action in gradients of Q, to produce a nonlinear 
sigma model. It turns out that performing this expansion requires a certain 
amount of care - the mathematical subtlety involved is known by the name 
of chiral anomaly. To do the expansion correctly, we have to resort to a vari- 



ant of nonabelian bosonization, which is the celebrated statement [14] that 
the free theory of n Majorana fermions has an equivalent representation by a 
level-one 0(n) Wess-Zumino-Novikov-Witten (WZW) model. Supersymme- 
try extends this to an equivalence between 2n + 2n Majorana fields (fermions 
and bosonic ghosts), and a level-one WZW model of fields taking values in a 
Riemannian symmetric superspace of type C\D, based on Mb x -^f where 
Mb = Sp(2n, C)/Sp(2n) and A^f = 0(2n). 

Using a supersymmetric extension of the bosonization rules for the chi- 
ral densities ip±ip^:, we are then able to compute the gradient expansion of 
the effective action. The Lagrangian of the CI | Dili nonlinear sigma model 
thus obtained contains a topological term, or winding number term, with 
topological coupling 6. Such a term is permitted by symmetry, since the 
massless Dirac Hamiltonian, Hq, depends on the choice of orientation of M^ 



or, in other words, reversing the orientation by a parity transformation, say 
by exchanging the two coordinates Xi and X2 of M^, takes Hq into an in- 
equivalent Hamiltonian. (Note that, quite generally, Dirac operators exist 
on spin manifolds, which are manifolds carrying a spin structure and thus 
an orientation.) The topological term is trivial for n = 1, but nontrivial for 
n > 2, since 

for n = 1 , 



n,(0(2„)/U(„))-, 2 to. „>2 

By reduction of the multi- valued action of the C\D WZW model, we show 
that the topological angle has the value 6 = Nir. 

The passage from random-mass Dirac fermions to the nonlinear sigma 
model is under good control, i.e. the terms omitted are small, for A^ ^ 1. 
The limit of a large number N of species is of course unrealistic. Fortunately, 
we can also control the case A^ = 2, if the kinetic energy is anisotropic, 
with the first species being more mobile in the xi direction and the second 
more mobile in the X2 direction. A closely related situation is relevant for 
the application of our results to disordered d-wave superconductors with a 
subcritical concentration of localized impurity spins ||T5|. For A^ ^ 1, or A^ = 
2 with large anisotropy, the nonlinear sigma model is at weak coupling, and 
the one-loop beta function predicts renormalization group flow to a Gaussian 
fixed point describing a perfect metal. The local density of states at E = 
in this case diverges logarithmically in the thermodynamic limit. We also 
compute the density of states for a finite system in the ergodic regime, and 
obtain agreement with the random- matrix prediction of 0. 

In contrast, for A^ = 1 the mapping on the nonlinear sigma model is far 
from being controlled. Even if we trust the mapping, the model is strongly 
coupled, and no safe statements can be easily made from it. We have to 
concede that our method fails in that case, and the problem is better analysed 
through its original formulation in terms of Dirac fields. Our conclusion 
thus is the one commonly accepted: random mass is a marginally irrelevant 
perturbation, and the theory in the infrared flows to the free-fermion point, 
where the density of states at zero energy vanishes. 

There now seems to be a conflict, but only superficially so, with Ziegler's 
rigorous proof of a positive lower bound for the density of states. The ap- 
parent discrepancy is resolved by observing that Ziegler works with a naive 
lattice discretization of the Dirac operator, thereby imposing an additional 
lattice symmetry on his model, as a result of which the Hilbert space de- 



composes into two decoupled sectors. In each of these, the Hamiltonian can 
be shown to be a pi-flux model (details of the argument will be published in 
a separate comment), which belongs to the time- reversal invariant Wigner- 
Dyson class Al and, in fact, is known [|1^, |1^ to have a finite density of states 
at zero energy. Hence Ziegler's rigorous result, albeit correct, is a statement 
about a different symmetry class, and does not falsify our results. We may 
safely ignore his proof in all that follows. 

There is another point of possible contention that deserves to be made 
clear. Recall that our large- A^ saddle-point analysis yields a target manifold 
consisting of two connected components. Thus there is a Z2 degree of free- 
dom in the effective theory, and we anticipate the existence of Z2 domains 
and domain walls, across which the nonlinear sigma model field jumps from 
one connected component of the target manifold to the other. (Note that 
in our computation of the one- loop beta function of the CI | Dili nonlinear 
sigma model, the nonperturbative effect of these domain walls is neglected.) 
The same Z2 degree of freedom arises for small A^, when the method of non- 
abelian bosonization is used. Indeed, nonabelian bosonization according to 



Witten |T^ transforms n Majorana fermions into a WZW model over the full 
orthogonal group 0(n), not just its identity component SO(n). The necessity 
to work with 0(n), which consists of two connected components, is readily 
seen by looking at the special case n = 1: a single Majorana fermion is not 
an empty theory, as would be the case with the trivial group S0(1), but is 
equivalent to a theory of local 0(1) = Z2 degrees of freedom, namely the 
2d Ising model. In a sense, then, passing from random-mass fermions to the 
effective description by a nonlinear sigma model, reintroduces the local Ising 
degrees of freedom underlying the free fermion. 

Having established the general formalism and its validity, we turn to a 
specific application: disordered non-chiral d-wave superconductors with bro- 
ken time-reversal and spin-rotation symmetry. Following a standard proce- 



dure [0, |T^, ^, we linearize the dispersion relation around four low-energy 
points (or "nodes") in wave vector space, and then add generic disorder (in 
the form of a random scalar potential, random complex order parameter, and 
random spin-orbit scattering) to place the model in class D. The number of 
species per node is A^ = 2 due to spin. By drawing on our previous results, 
we map the low-energy physics of the quasiparticles on a nonlinear sigma 
model. Doing this for every node separately, we obtain kinetic energy terms 
that are spatially anisotropic, and a topological angle which is positive for 
two of the nodes and negative for the other two. The inclusion of quasipar- 

8 



tide scattering between the nodes merges the four anisotropic models into a 

single isotropic theory. In the last step, all the topological terms cancel, as is 

required for a non-chiral superconductor invariant under reflections of space. 

Next, we recall the 2d phase diagram proposed by Senthil and Fisher 



111 , and address the puzzle why the two stable phases of the 2d random- 
bond Ising model seem to have three counterparts in class D. To resolve 
this mismatch, one needs to understand the nature of the multicritical point 
A1* and the renormalization group flow in its vicinity. Our resolution of the 
puzzle is very simple: the "multicritical" point A1*, where the three phases 
of class D meet, is to be identified with the free-fermion point! The evidence 
in favor of such an identification is twofold. Firstly, by an argument involving 
nonabelian bosonization, we show that the free-fermion point sits right on 
the boundary of the 2d metal phase. Secondly, we argue that (barring vortex 
disorder) the thermal metal phase of class D can only be realized in the 
enlarged parameter space available for more than one species (A^ > 1), and 
is absent in the fundamental case A^ = 1. This makes it natural to invert 



some of the flows assumed in []rT| and arrive at a plausible 2d phase diagram, 
where the paramagnetic and ferromagnetic states of the Ising model match 
the thermal insulator and thermal quantum Hall fluid phases of class D for 
A^ = 1, and there is no redundant third phase. 

Another suggestion by Senthil and Fisher ITTl says that, from the per- 



spective of the nonlinear sigma model, the existence of the insulating phases, 
which are distinguished by a quantized thermal Hall conductance, can be 
attributed to the presence of the topological term in the action functional. 
One might then think that this very term is what is driving the phase tran- 
sition from the thermal insulator to the thermal quantum Hall fluid. That 
scenario is widely accepted for a close cousin, namely Pruisken's nonlinear 
sigma model [^] of the integer quantum Hall effect. However, as was stated 
earlier, the topological term of the CI | Dili nonlinear sigma model is trivial 
for n = 1. Since the phase transition occurs for all n, including n = 1, it is 
hard to see how the topological term could be its driving agent. A better 
scenario is to attribute the phase transition to the Z2 degree of freedom of 
the CI I Dili nonlinear sigma model. 

Finally, we delve into the quasiparticle physics of chiral p-wave supercon- 
ductors, to shed more light on the effect of vortex disorder in those systems. 



Read and Green ||T3l have argued, both intuitively and formally, that ran- 
domly placed vortices are a relevant perturbation at the free-fermion point 
of a superconductor with (mean-fleld) order parameter symmetry Px+iy Un- 



like the perturbation by a random mass matrix, which can become relevant 
only in the extended parameter space available for N > 2, vortex disorder 
turns out to be relevant even in the fundamental case of spinless particles 
(A^ = 1). Read and Green further propose that vortices drive the supercon- 
ductor from the free-fermion point into a thermal metal phase, the effective 
field theory for which lacks, as it stands, the Z2 (or Ising spin) degree of 
freedom of the nonlinear sigma model we obtain for random-mass fermions. 
We believe that their proposal is correct, and in order to provide supporting 
evidence for it, we outline an indirect argument, passing through a variant of 
the Chalker-Coddington network model with local 0(1) invariance. (We are 
unable at present to give a direct field-theoretic proof, as we do not know how 
to incorporate vortex disorder into the nonabelian bosonization scheme.) 

What we are learning, then, is that vortex disorder exerts a drastic in- 
fluence on the effective field theory: while there exists a local Z2 degree of 
freedom (and domain walls in it) when vortices are absent, this degree of 
freedom is suppressed when vortices are inserted. Thus the cases with and 
without vortices map on different field theories. Actually, there must exist 
a whole one-parameter family of such theories, as the suppression proceeds 
continuously. Motivated by random-matrix limits, we propose to refer by the 
(hybrid) name BD to the generic symmetry class including a variable amount 
of vortex disorder. Class D is then viewed as a subclass of the generic class 
BD. Of course, the field-theoretic distinction between classes does not fail 
to leave its imprint on the quasiparticle physics. We expect that when all 
quasiparticle states are localized, i.e. in the thermal insulator and quantum 
Hall fluid phases, the local density of states at zero energy vanishes in the 
absence of vortices (class D), but becomes flnite when an extensive number 
of vortices is inserted (class BD) . 

In contrast, the physical distinction between systems with and without 
vortices is of relatively minor consequence in the metallic regime of extended 
states. There, and for the case of a flnite system and energies E below the 
Thouless energy E^h, the nonlinear sigma model can be evaluated in zero- 
mode approximation. For the total density of states in the two cases at hand, 
we obtain: 

Pd{E) =u+ ^^^ ^ , Pbd{E) = u+ 16{E) , 

where the scale of variation is set by u, the inverse of the level spacing 
for energies much greater than the mean level spacing (but still less than 
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i?Th) • The first expression results on integrating over the full target manifold 
(with Z2 present), the second from the reduced target (with Z2 completely 
suppressed). These results coincide and, as we shall see, not accidentally so, 
with the large- A^ limits of the density of states for the Haar random-matrix 
ensembles on S0(2iV) and 0(2A^), respectively. 

The material of the paper is arranged as follows. In Sections ^^, we 
introduce a Hamiltonian for N species of Dirac fermions with random mass, 
and set up a field-theoretic representation of its Green functions by a su- 
persymmetric Gaussian functional integral. Renormalization of the disorder- 
averaged field theory generates a total of four marginal perturbations, the 
flow equations for which are worked out in Section ^ Guided by this, we set- 
tle for a specific choice of Lagrangian and proceed to map it on the CI | -Dili 
nonlinear sigma model, in the limit A^ ^ 1. The composite field Q is in- 
troduced and subjected to a saddle-point approximation in Section |, while 
the structure of the saddle-point manifold is elucidated in Section ||. Section 
^ is the longest of the paper. There, we expand the effective action for Q 
in gradients, using a supersymmetric extension of the method of nonabelian 
bosonization, which is developed in two major subsections. In Section H, we 
renormalize the nonlinear sigma model (neglecting the Z2 degree of freedom) 
and show that its coupling flows to zero, resulting in a perfect metal with 
infinite thermal conductivity. We compute the density of states for a finite 
system, as well as in the thermodynamic limit. In Section ^ we clarify the 
relation between the model considered and disordered d-wave superconduc- 
tors. Section ^ assembles various arguments that focus on the role of the 
free-fermion point inside class D, and culminate in the proposal of a local 
phase diagram. Finally, in Section [T^ we comment on the role of vortices 
and the significant modification they cause in the nonlinear sigma model. 

To avoid confusion, we emphasize that the present paper will always be 
concerned with the physics of Majorana fermions. Nevertheless, the words 
"Dirac fermions" or "Dirac theory" will be frequently encountered, the reason 
being that we carry out the disorder average by using the supersymmetry 
method, and a single Majorana fermion does not have a bosonic analog. 
The standard trick to get around this problem is to "square the partition 
function" , making two Majorana fermions out of one, and then combine the 
two to form a Dirac fermion. The latter can be augmented by a bosonic 
b-c ghost system to make the theory supersymmetric and cancel all vacuum 
graphs. 
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2 Dirac Hamiltonian with random mass 

The first-quantized Hamiltonian H for A^ species of Dirac fermions with 
random mass in two dimensions is written as 

^ ~ \-2id -M^{x)) ■ ^"^^ 

Here M = {Mki) is an N x N mass matrix, and we use the convention 
d = |(i9i — id2), B = \{di -\- 182) ■ Aside from being Hermitian {W = H), 
this Hamiltonian enjoys an important symmetry property: 

H = -aiH^^ai , 

which will be called particle-hole symmetry, or the symmetry of class D. 
The superscript T denotes joint transposition in particle-hole and position 
space (note d"^ = —d). The particle-hole symmetry of H is dictated by the 
Lagrangian for Majorana fermions, 

(Summation over repeated indices is understood.) An immediate conse- 
quence of the particle-hole symmetry is that the eigenvalues of H occur 
in pairs with opposite sign ±E. Hence, the density of states is symmetric 
with respect to the point E = 0. By the same token, there exists a relation 
between retarded and advanced Green functions G'^{E) = (i? ± iO — H)~^: 

G-{E) = -aiG+{-Efai . 

When the mass matrix is a constant multiple of unity, M{x) = m x Ijv, one 
easily computes the local density of states: 

u{E) = {iiTeii)-\-HmTT{E-iO-Hy^ 

2N ^ ^ [ <Pk E-ie 

lim Im 



TT e^o+ J (27r)2 {E - it)^ - k^ - m"^ 



N\E\,,^2 2 



2n 
which vanishes linearly at i? = in the limit of zero mass. 
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In what follows, we take the entries of the mass matrix to be Gaussian 
distributed random variables with zero mean and second moments given by 



{M,,{x)Mki{y)) = {2gM/N) 5u 5,k 5(x - y) . 



(4) 



The density of states on average over the fluctuating mass matrix can be 
obtained by computing the average Green function. This will be done by an 
adaptation of Efetov's supersymmetry method ||22[. Although that method 
is in principle quite standard, we will see that its application to the present 
situation features some peculiarities which are well worth explaining. 



3 Supersymmetric integral representation 

We are going to employ a supersymmetric integral representation to set up 
the calculation of the disorder averaged Green functions. Although we con- 
centrate on the case of a single Green function {n = 1) for notational simplic- 
ity, the generalization to arbitrary n> 1 will be immediate |H] . We introduce 
2N supermultiplets (pk and (^^ (fc, Z = 1, ..., N), each containing two fermionic 
components ■ip± and two bosonic superpartners h and c.0 Each supermultiplet 
is arranged in the form 



c 
V b ) 



c 



where the index / 
of is defined by 



A^ was omitted. The orthosymplectic transpose 



0* = (7/.+ ,V^_,6,-c), (5) 

and (f)^ (defined similarly) is the orthosymplectic transpose of cj). Note that 
this choice of transpose determines an inner product which is skew: 



(p^cj) = ip+ip- + 4'-4'+ + be — cb 



(6) 



Moreover, the orthosymplectic transpose on defines a compatible transpose 
of 4 X 4 supermatrices T by the requirement (T0)* = 0*T*. 

''We use the conformal field theory convention of denoting antiholomorphic fields by 
an overbar. To avoid confusion where it might otherwise arise, complex conjugation is 
denoted by an asterisk * instead of the overbar. 
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The field-theory Lagrangian for A^ species is taken to be 

Lo = 0190, + 0190, + 20lM,,0, - zi?0lS30, , (7) 

where 

S3 = Isusy ® 0-3 = diag(+l, -1, +1, -1) , 

and summation over repeated k, I indices is understood. If j and j are source 
fields, the Green functions of the Dirac Hamiltonian (^ are generated by the 
functional 

Z[j] = ( y P0 1)0 exp - y d^x (Lo + 0V + 0*j) 

as is easily verified from the fact that Lq is the quadratic form constructed 
by sandwiching the operator 

■rrr m^ , f iM - lE 2(9 \ 

1,„,, ® t{H -E) = l,,3y ®y 2d -zM^ -ze) 

between 

{ij+,b,ip+,b) and y 

V -cj 

and integrating by parts to symmetrize the terms with derivatives. In order 
for the Gaussian functional integral over 0, to make sense, the bosonic fields 
must be related to each other by complex conjugation *: 

&; = q, c; = 6, . (8) 

Given this convention, the functional integral is purely oscillatory for E E'R, 
and converges for energies in the upper half plane (Imii^ > 0). As stated in 
the introduction, the Lagrangian Lq in the limit E ^ acquires an invariance 
under global transformations, 

0,(a;) f->T-0;(x) , 0/(x) H^T-0,(x) , 
0Kx) H. 0Kx) ■ T\ 0Kx) ^ 0Kx) ■ T* , 

for elements T G OSp(2|2), i.e. if T obeys T^T = 1, with T* being the 
orthosymplectic supermatrix transpose defined above. 
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On performing the disorder average specified by (Q), the Lagrangian for 
the case E = becomes 

Li = 0190, + <p\d<P^ + (gM/N)^^'^ , $(1) = fy,4>'<l>k ■ (9) 

By simple power counting, the operator $*^^^ is a perturbation of the massless 
Dirac theory which is marginal in the renormalization group sense. ^^^^ may 
be marginally relevant or marginally irrelevant. The question of which is the 
case, is decided by calculating the short-distance expansion of the operator 
product of $''^-' with itself. 



4 Marginal perturbations 

It turns out that the operator product expansion of ^^^' with itself generates 
three additional operators: 



$(2) 
$(3) 
$(4) 



A(t>k(t>\(t>i , 



for A^ > 1. Each of these is marginal and invariant under OSp(2|2), and 
the set $(^), ..., $*^^^ exhausts the set of marginal perturbations permitted 
by that symmetry. Note that the situation simplifies for A^ = 1, where 

$(1) = $(3) ^ $(4)^ g^n^ $(2) ^ Q f^Qj^ (^|^_ 

Given that the operators $^"'' (a = 2,3,4) are generated by the renor- 
malization group flow anyway, the natural procedure is to include them in 
the bare theory. Our goal thus is to renormalize the extended Lagrangian 

4 

L2 = ^Id^i + <Pld<Pi + iV-1 J2 9a<^^"^ ■ (10) 

To do so, we use the fact p3| , p4| that, if the leading short- distance singularity 



of the operator product expansion (OPE) is 



the one-loop beta functions determining the RG flow with increasing cutoff 
length scale i, are given by 

2 . _ dg^ 7T 



' ■ "'9-y "^ \~^ r^a0 
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The expansion of the operator product $*^")(a;)$^^-'(x') is completely deter- 
mined by the OPEs for holomorphic {z = Xi + 1X2) and antiholomorphic 
{z* = Xi — 1x2) fields. These are 

1 r —ah _ _ 1 A T-"^ 

rkizWw) ~ -^ ^^^^ , 0^(z*)0fK) ~ _L^^[^ , (11) 



where 

r = o"! (g) Eff + io-2 ® Ebb 



/O 1 0\ 

10 

1 

\0 -1 0/ 



determines the orthosymplectic structure. Using these formulas, a lengthy 
but straightforward calculation yields 

9i = 9i92 - N'^ {.9192 + 9i9i - 9293 + 9^94.) , 

92 = Udi + 92) + N'^ {9193 + 9i94 - 9I - 9394) , 

93 = -{91+ 92)93- N'^ {9I- 9i92 + 9i94 + 9293 + 9I + 9394) , 

94 = -N'^ {gig2 + s-ifi-a + 9293 + 9T} ■ (12) 

For A^ = 1 these equations simplify greatly. In that case, $''^-' = 0, so 
the coupling g2 must be dropped, and since ^^^' = ^^^^ = ^^^\ the only 
coupling to renormalize is Qm = 9i + 93 + 94- By linearly combining the 
above equations, we get 

dgM _ _'2gli f-,o\ 

dlni~ TT ' ^^"^^ 

which reproduces the well-known result []I| that random mass is a marginally 
irrelevant perturbation (for A^ = 1). 

Let us quickly review how to deduce from this result the behavior of the 
local density of states v^E^qm) at small energies E and weak coupling qm- 
We start from the formula 

u{E,gM) = T^-'Re{{ilJ-^+ + ^-^+m) , 

where the functional average is defined w.r.t. the Lagrangian (^ for A^ = 1. 
The OPE between $*^^^ and the operator that couples to E^ $*^'') = 0*^30, 
reads $(i)(a;)$(°)(0) = -(27r2)-i$(o)(0)/|a;|2 + .... It then follows, by the 
same principle we just used for the coupling qm-, that the operator $*^°) has 
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scaling dimension 'J^Qm) = 1 — gu/T^ + 0{.g\i)i and the RG flow equation for 
E is dE/dlni = (2 — 'y{gM))E. By making in the functional integral an RG 
transformation which changes the cutoff from £o to £, we obtain 



We know from ([T3| ) that the coupling gu^d-) flows to zero with increasing L 
We also know the density of states of the pure system to be independent of 
the cutoff: u{E,Q;t) = \E\/2tt. Hence, by integrating the flow equations for 
gM and E up to the infrared cutoff, which for small enough energies is given 
by the system size L, we obtain 

iy{E,gM;io) = P(l + (2W^)ln(V4)) • 
/vr 

We see that, for A^ = 1, the density of states at ii^ = remains zero, and 
randomness in the mass only changes the slope of t/ ~ \E\ by a logarithm. 

For A^ > 1 the flow equations for the couplings gi, ...,g4 are not very 
transparent, and the precise form of the RG flow remains unclear in general. 
However, for A^ ^ 1 the equations reduce to 

gi ± g2 = ±l{gi ± g2)'^ , 93 = -{91 + 92)93 , 94 = 0. 

Clearly, the parameters g^ and g^, if initially set to zero, remain zero under 
the flow. The relevant couplings are gi and g2, which are seen to be attracted 
to the line g = gi = g2, with g satisfying the equation 

"' m = +^ . (14) 



d\ni TT 

We observe that the sign on the right-hand side has been reversed as com- 
pared to A^ = 1, and the coupling g is now relevant. 

These considerations motivate us to adopt the modified Lagrangian 

L3 = ^'Mi + <P\d<P, + {g/N) ($(1) + $(^)) . (15) 

At the level of the Dirac Hamiltonian @, the inclusion of the operator $*^^) 
corresponds to adding a term 



H ^ H 



A 
At 
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where A is a skew-symmetric N x N matrix whose entries are Gaussian 
distributed random variables with zero mean and second moments 

{Ai,{x)AM{y)) = {2g/N){6ikSji - 6ii6jk) Six - y) . 

Such a term is permitted by the fundamental particle-hole symmetry (|^) of 
class D. Note that for A^ = 1, skew-symmetry forces A to vanish identi- 
cally. ( Equivalent ly, $''^'* =0 for A^ = 1, and L3 coincides with the original 
Lagrangian.) For the application to disordered superconductors (Section ^, 
we put N = 2. The indices i,j, k, I then label the spin degrees of freedom of 
the electron, while the random functions M{x) and A(x) acquire a physical 
meaning as random scalar potential, random spin-orbit scattering, and the 
random part of the superconducting order parameter (or, to be precise, linear 
combinations thereof). 

Before embarking on our project of analysing L3, we wish to mention 
three other special choices of the couplings, where we can easily understand 
the nature of the RG flow. The first case is gi = gs = 0, N > 2. By 
tracing the couplings back to a Hermitian random Hamiltonian, we find the 
constraint (71 ± (74 > 0. Thus, setting gi = forces g^ = 0. Inspection shows 
that the random Hamiltonian in this limit has a higher degree of symmetry: 

H = -(ai ® lN)H^iai (g) l^v) 
= +((T2 ® lN)H^{(r2 ® In) , 

where the second equality is a kind of time-reversal invariance, placing H in 
class Dili [0. The only coupling that remains is g2, which flows according 
to g2 (X {N — 2)gl and thus is relevant for N > 2. The methods developed 
in the present paper are readily adapted to deal with this case. 

The second special case is obtained by setting A^ = 2 in the previous one. 
The underlying random Hamiltonian then is of the form 

TT _ ( —2id + ao2 

~ \ -2id + aa2 

Because this commutes with diag(o"2, cr2), the Hilbert space decomposes into 
two sectors invariant under the action of H. (For N > 2, the Hilbert space 
does not decompose.) The Hamiltonian restricted to either sector is a Dirac 
fermion in a random abelian vector potential ±a{x), which is class Alll [j^ ^, 
and explains p5l why the single coupling g2 is exactly marginal for N = 2. 



The extra time-reversal symmetry here does no more than connect the two 
restrictions oi H hj a Lie algebra homomorphism. Apart from causing every 
energy eigenvalue to be doubly degenerate, this is of no consequence, and the 
symmetry class is AIII instead of Dill. 

There exists a third special case where the flow equations close on a single 
coupling. This is g2 = g-i = = gi + g^ with gi > 0, and again N = 2. Here 
the underlying random Hamiltonian is of the form 

rr^ f V(^2 -2td 

\ -2id Va2 

By a simple conjugation H' = UHU^^, this can be transformed into H' 
diag(if+, H_) where 

±V -2id^ 



^^ ' -2id ±v 

are two copies of a Dirac fermion with random scalar potential V{x). In the 
transformed basis, each copy satisfies 

H± = a2Hla2 , 

which is the defining equation of the "symplectic" Wigner-Dyson class All. 
The extra particle-hole symmetry (class D) here does no more than relate the 
first copy to the second copy by a Lie algebra homomorphism. It is therefore 
redundant, and the symmetry class remains All. 

By numerically solving the flow equations (|12D, we find that the line 
92 = 93 = = gi + di has a large basin of attraction for N = 2. In particular, 
it attracts the initial condition gi = g2, g^ = g^ = in that case. 

According to [^, the classes D, Dill and All exhaust the set of symmetry 
classes whose field-theory representation is acted upon by the present form of 
OSp(2n|2?T,) (preserving a symplectic structure for bosons and an orthogonal 
structure for fermions). We therefore expect no further accidental closures 
to occur. 



5 Composite field and diagonal saddles 

As we saw, for A^ ^ 1 the coupling g = gi = g2 is certain to increase 
under renormalization, leading to a strong-coupling problem in the infrared. 
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To handle this problem, it is mandatory to switch to a dual formulation in 
terms of a Hubbard-Stratonovich field Q ~ 0;0f + 0^0;. 

The composite field Q is brought on stage by doing two Gaussian integrals 
in a row. To prepare these steps, we reorganize the Lagrangian (^) as follows: 

L, = 0190, + 0150, + ^ STr (0,0*, + 0,0*,) (0,0f + 0,01) . 

The quartic interaction can now be decoupled in a first step by introducing 
an auxiliary 4x4 supermatrix field Q: 

U = 0190, + 0190, + STrg(0,0l + 0,01)-^ STr g2 



**^»(n)a;)-^--- 



STrg2. (16) 



In the second step, one does the Gaussian integral over the Dirac field, re- 
sulting in the action functional 

S4Q] = -^J'f^STrQ-' + f STVlu (^f ^^^) . (17) 

where STr combines the operations of taking the supertrace STr and inte- 
grating over position space. The factors E3 under the logarithm appear after 
the transformation 0, \-^ ^30, and 0* ^-^ 0*^3, correcting for the fact that 
is not the complex conjugate of 0*; see the definition (|^) and the relations 
(H). The above result is for E = 0. To restore the dependence on energy, we 
need to shift Q ^ Q — {iE/ 2)113 under the argument of the logarithm. Note 
that all manipulations done so far were exact. 

As a result of the transformation to Q, the parameter A^ now appears as a 
factor in the exponent, while Q itself is ignorant of the number of species. For 
large A^, this suggests treating the Q integral in saddle-point approximation, 
which is what we are going to do next. (Note, however, that the saddle-point 
approximation is uncontrolled for A^ = 1.) The saddle-point equation reads 

Q{x)/g = l{x\{Q- dQ-^d)-^\x) + {d ^ B) . 

We look for a spatially homogeneous solution of the form Q{x) = fiT,^. The 
saddle-point equation then reduces to 
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Cutting off the integral in the ultraviolet by \k\ < I/Cq yields the equation 

7r/(7 = ln(l + (2/i£o)"') , 

and by inversion, 

/x = (2£o)-VVe-/.-l. (18) 

For weak disorder {g <C 1), this is well approximated by /x ~ (2£o)~^e~''/^^. 
From formula (|1^ we then infer that /i obeys the renormalization group 
equation 

"^ +/3(^)A^^(£,^) = 0, 



^ (i In £ dg ^ 

and thus has the meaning of a dynamically generated mass. In other words, 
the mass scale /x remains invariant: 

M4,^(4)) = M^,^W), (19) 

under an RG transformation Iq v-^ L This means that the Dirac field under- 
goes the phenomenon of dimensional transmutation: while the bilinear 0^0; 
has dimension one at the free-fermion point, its dimension approaches zero 
on large scales. At the same time, Q ~ (f)i(f)\ + (f)i(f)\ acquires a nonvanishing 
expectation value. 

In principle, all diagonal matrices Q with entries ±/i are candidate so- 
lutions of the saddle-point equation. It may happen, however, that some of 
these do not lie on the integration domain, or cannot be reached by analytic 
continuation from that domain. Therefore, in order to be able to select the 
proper solutions, we have to be more explicit about which integration domain 
to choose for the superfields 0;, 0/ and Q. 

6 Saddle-point manifold 

The issue of how to choose correctly the integration domain for the superfields 
0,0 and for the supermatrix Q was carefully addressed in |^. Drawing on 
(but not repeating) the detailed analysis of that reference, we argue as follows. 
Let Imii^ > for definiteness. Then, in order to have convergent integrals 
over the bosonic ghosts of the Dirac version of the theory, we must impose 
the conditions (H) or, equivalently, 4>aB — 4>\B^'i- Note that there is no need 
to require such a relation among the Grassmann fields, and we do not impose 
any such requirement. 

21 



Next one easily verifies that the bihnear A = (f)i(f)\ + (f)i(f)\ in the Lagrangian 
L4 is odd under the orthosymplectic transpose: A^ = —A, so A G osp(2|2), 
by definition of the orthosymplectic Lie algebra. Via its coupling to A in L4, 
the supermatrix Q inherits the same property: 

Q = -Q' e osp(2|2) . 

Now, in the process of decoupling the interaction $*^^^ + $*^^\ severe con- 
vergence problems arise in the boson-boson (BB) block denoted by Qbb- 
However, close inspection shows that the integrals over Qbb can be arranged 
to exist (without ruining the convergence of the integrals over the bosonic 
Dirac ghosts b, c) by choosing the following parametrization: 

gBB = V' + 7e''S3e-^, 

where X, F G osp(2|2) are odd resp. even with respect to conjugation by S3: 

S3XS3 = —X , S3FS3 = Y , 

and the even (or bosonic) parts Xq, Yq of X, Y are subject to 

The parameter 7 > is in principle arbitrary but, anticipating its saddle- 
point value, we set it to 7 = fi. This deals with the BB sector. The situation 
in the FF sector is more benign. There, convergence of the integrals over 
Qpp is simply achieved by requiring Qff = Qff- 

Now we can compare the above parametrization for Q to the diagonal 
solutions of the saddle-point equation. The structure of the solution in the 
BB block is uniquely fixed by the convergence requirements to be 

Qo,BB = /^0"3 • 

Other diagonal solutions cannot be reached by deformation of the integration 
manifold without crossing some singularities of the integrand. 

In contrast, in the FF sector analyticity provides no selection rule on the 
saddle points; there, the integrand does not possess any poles, but only has 
zeroes as a function of Qff, so that the path of integration can be analytically 
continued to cross any saddle point in that sector. However, in the weak- 
coupling (or large A^) limit, some of the saddles contribute in a negligible 
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manner. The dominant contributions come from those that minimize the 
super- dimension of the transverse manifold. This extremahty condition is 
fulfilled when the positive and negative eigenvalues of Qcff are equal in 
number. In the present case, this criterion leaves two possibilities: 

<5o,FF = ±/^o-3 • 
In summary, we retain as dominant (diagonal) saddle points: 

Qo = l^ (±^FF + ^bb) ® (X3 . (20) 

Recall now the existence of a global G = OSp(2|2) symmetry (at E = 0). 
The symmetry group acts on the supermatrix Q by conjugation: Q{x) i— ^ 
TQ{x)T~^ {T ^ G). Clearly, such transformations leave the action functional 
S^IQ] in (^) invariant. As a result, the saddles of S^IQ] are degenerate: given 
any solution Qq of the saddle-point equation, we get an orbit of solutions by 
acting on Qq with the symmetry group G. Because the stability condition 
hQoh"^ = Qq divides out a subgroup H = GL(1|1), called the stabilizer of 
Qo, the orbit of G on Qq is a coset space G/H = OSp(2|2)/GL(l|l). Notice 
that this is the orbit of the complex symmetry group. 

To carry out the saddle-point integral, we need to restrict the bosonic 
degrees of freedom of G/H to a real submanifold. We now briefly describe 
this submanifold, setting the Grassmann variables temporarily to zero. The 
bosonic part of the complex supergroup OSp(2|2) is 0(2, C) x Sp(2,C). In 
the BB sector, the symmetry group is Sp(2, C), which arises as the group of 
transformations T, 

leaving invariant the symplectic form 

1\ fb 



be- cb ,>.,^, , , „ , . 

1 Oy \c 

Those symplectic transformations of (6, c) and (6, c) that preserve the re- 
ality conditions (||), form a real subgroup Sp(2,M). Division by the sta- 
bilizer, which is isomorphic to U(l) G H, then yields the BB manifold 
Mb = Sp(2,M)/U(l) c^ H^. This is a two-hyperboloid. (For general n, 
we get Sp(2ra, M)/U(n), which is a noncompact symmetric space of type CI.) 
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In the FF sector, the symmetry group acts by 0(2, C), which is under- 
stood to be the invariance group of the symmetric form (or "orthogonal" 
structure) 

V'+^_ + ^_^+ = (V'+,^-)(j o){t~ 

We have emphasized the crucial fact that we never impose any reality condi- 
tions on the fermions. What determines then the real FF subgroup? The an- 
swer is that reality here enters through the Hermiticity constraint Qff = Qff 
(needed for convergence), which selects from the complex symmetry group 
0(2, C) a real subgroup 0(2), the usual orthogonal group in two dimensions. 
Dividing by the stabilizer U(l), we obtain Mp = 0(2)/U(l) ^ Z2. (For gen- 
eral n, we get 0{2n)/U{n), a compact symmetric space of type Dill.) Thus 
the base of the saddle-point manifold is M-q x Mp = H^ x Z2. On reinstat- 
ing the Grassmann variables, we arrive at a saddle-point supermanifold, Xi, 
which is a Riemannian symmetric superspace of type CI | Dill 0. 

A distinctive feature of Xi is that, instead of being connected, it consists 
of two disjoint pieces. This fact was already discovered in |p, and we are 
now going to elaborate briefly. The special orthogonal group SO (2) acts on 
two Majorana fermions (^, rj) as 

^\ ^ /cos9 -sm9\ n 
7] J y sin 9 cos 9 J yi] 

while the corresponding SO (2) action on the Dirac fermion ip± = ^ ± ir^ is 

^+\ fe^' \ /^H 

This does not exhaust the symmetries of the FF sector. The Majorana 
fermion is known to possess an additional discrete symmetry, which is reflec- 
tion ^ *-^ rj (an element of 0(2) with determinant minus one), or in Dirac 
language: 

V,J"(-. o)[i.J=-'"'[r 

Consider now the action of these symmetry transformations on the saddle 
point Qo.FF = CTs in the FF sector. The action of SO (2) simply fixes the 
saddle point: e^^'^^ase"*^""^ = era, and thus is trivial, while the discrete 0(2) 
transformation by (72 reverses the sign: 0-3 \-^ cr2cr30"2 = —(^s- The latter then 
explains the origin of Z2 in the symmetric superspace Xi. 
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In summary, all of the saddle points of 5*5 [Q] lie on a single orbit of the 
bosonic symmetry group 0(2) x Sp(2,R) C OSp(2|2), but the orbit is not 
connected, as 0(2) comes in two pieces. We mention in passing that this 
peculiar feature of universality class D also occurs in class Dili p . 

Although we have focused on the formalism for a single Green function 
[n = 1), all of our considerations readily extend to the case n > 1. Supervec- 
tors and supermatrices simply become longer or bigger. The symmetry group 
inflates to G = OSp(2n|2n), the matrix E3 tensors to S3 = Isusy ® o"3 1„, 
and the stabilizer becomes H = GL(n|n). The order parameter space for a 
general value of n is denoted by X„, and has bosonic submanifold 

Mb X Mf = (Sp(2n,M)/U(n)) x (0(2n)/U(n)) . 

Again this consists of two disjoint pieces, corresponding to the two connected 
components of 0(2n) (with determinant plus or minus one). In the following 
section, we shall work with the extension to general n. 

As an important corollary to the discovery of a Z2 degree of freedom in the 
saddle-point manifold, we anticipate the existence of Z2 domains and domain 
walls, across which Q jumps from one connected component of the saddle- 
point manifold to the other. We expect this Ising-like degree of freedom to 
be very important for the phenomenology of the insulating phases of class 
D (Section |1^, where Q fluctuates strongly. On the other hand, when the 
field is stiff, domain walls are costly in energy, and therefore they should be 
of minor relevance in the metallic limit we shall focus on in Sections 0-|^. 

7 Gradient expansion 

To summarize the current state of affairs, the problem at hand has a global 
G = OSp(2?2|2n) symmetry, which is broken to GL(n|n) by making a naive 
saddle-point approximation. Whether the symmetry is truly broken or not, 
will have to be decided at a later stage. (The Mermin- Wagner- Coleman the- 
orem, stating that continuous symmetries of compact type cannot be broken 
spontaneously in two dimensions, does not apply in the present case, as 
the saddle-point manifold is a noncompact superspace. Hence, the question 
whether symmetry breaking occurs or not remains open for now.) Our saddle- 
point analysis has identified an order parameter Q = fiq, where q = TJ^^T'^ 
takes values in a symmetric superspace X„. The low-energy configurations 
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of the action S^ [Q] in ([T7| ) are given by slowly varying fields 

q{x) = T(x)S3T(x)"^ (T(x) G G) . 

These are the Goldstone modes of the broken G symmetry. Note that q lies 
(on an adjoint orbit of G) in osp(2n|2?T,), and satisfies the constraint q^ = 1. 
Our next goal is to derive the low-energy effective action for the Goldstone 
modes q{x). (We shall neglect field fiuctuations transverse to the saddle-point 
manifold, since these are massive. We shall also assume q{x) to be smooth, 
which means we will ignore domain walls in the Z2 degree of freedom.) On 
general grounds, the effective action must be of the form 

Scs [q] = T^—F / d'^x STr df,q d^q + — / d^x e^„ STr q d^q d^q . (21) 

The two terms displayed are the only ones that contain no more than two 
derivatives, respect rotational invariance of position space, and are compati- 
ble with global G symmetry. (Later we will add a symmetry-breaking term 
for finite energy E' 7^ 0.) The second term is a topological or winding number 
term. It arises by pulling back the closed two-form STr g dg A dg (called the 
Kahler form) of X„, and is nontrivial since 

n2(X„,) = n2(MB X Mf) = TI2 (0(2n)/U(n)) = Hi (U(n)) = Z 

for n > 1. (In contrast, the n = 1 winding number n2(0(2)/U(l)) = 
is trivial. It may seem curious that there exists such a marked difference 
between n = 1 and n > 1, as the topological term is supposed |]lll to play an 



important role in controlling the phase transition between the two insulating 
phases of class D in two dimensions. We will suggest an explanation later, 
when we discuss the phase diagram.) Quantitatively put, for any closed 
surface S we have 



-— / (fxenuSTiqdnqduq E Z 
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For n = 2 the correct normalization factor (327rz)^^ can be figured out by 
direct calculation, using 0(4)/U(2) ~ S^ x Z2. The generalization to ra > 2 
is dictated by the connection with the multi- valued action r[M] in equation 



(p^) below. We note that the topological term is odd under parity xi <-> 0:2. 



Nevertheless, its presence is permitted, as a choice of pure Dirac Hamiltonian 
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Hq = cTipi + (J2P2, as compared to Hq = o"2Pi + CT1P2, implies a choice of 
orientation of position space. 

Our task now is to work out what the values of the couplings / and 6 
are. This is easily done for /, by straightforward gradient expansion of 

.„l4.„SDet(''t* ,,!)■ (^^) 

On the other hand, the value of the topological angle 6* is a subtle issue. 
This is not easily found by gradient expansion. The reason is that topological 
excitations such as instantons, whose topological charges are what is counted 
by the winding number term, necessarily involve large variations of the field q, 
thereby defying simple considerations based on Taylor expansion of the action 
functional. One viable option would be to evaluate both Ses[(l] and 5*5 [g] on a 
well-chosen configuration q{x), say an instanton solution (such solutions exist 
for n > 1) and equate the two answers to determine 6. Unfortunately, the 
evaluation of 5*5 [g] on an instanton is neither easy nor rewarding (which is to 
say that it does not lead to any significant payoff beyond solving the technical 
problem at hand). A powerful standard trick for computing determinants 
of Dirac operators is to take the logarithmic derivative with respect to a 
parameter s, and integrate over s at the end. Such a strategy is doomed to 
fail here, as differentiating a winding number necessarily produces zero (and, 
moreover, topologically distinct sectors cannot be continuously connected 
with one another). 

The most elegant and painless scheme for extracting the low-energy effec- 
tive action from S^ [q] proceeds via the method of nonabelian bosonization. 
Following a celebrated paper by Witten [H, this method has become a stan- 
dard tool for dealing with perturbations of the free fermion theory. Unfortu- 
nately, in the supersymmetric context of Dirac fermions augmented by a b-c 
ghost system, the method is less established. Therefore, in the subsections 
that are appended below, we will make a digression to explain the needed 
extension of nonabelian bosonization, using the language of functional inte- 
grals. For now, we describe how the method is used to achieve our goal of 
expanding 5*5 [g] in gradients. 

We start by undoing the integration over the fields (p, 0: 



e 



s^kW ^ T>(PV(P exp - / d^x (0*(90 + 0*(90 + /x0*g0 + /x0*g0) 



The notation is the same as before, except that we have dropped the sum- 
mation over the species index (/ = 1, ..., A^) and divided 5*5 by A^ accordingly. 

27 



The next step is to bosonize, using the principle of Bose-Fermi equivalence 
in two dimensions. Witten's nonabelian version asserts that the free theory 
of 2n Majorana or n Dirac fermions has an equivalent representation by a 
Wess-Zumino-Novikov-Witten (WZW) model with target 0{2n) resp. U(?t,) 



at level /c = 1. As will be shown in Subsections |7.1| and [7. 2| below, this equiva- 
lence can be generalized to allow for the presence of a h-c ghost system on the 
free-fermion side. The equivalent "bosonized" theory turns out to be what 
might loosely, but only loosely, be called an "OSp(2n|2n)" WZW model. We 
will give the precise definitions later. Here we simply state that the WZW 
field, which we denote by M , takes values in a subspace of the complex su- 
pergroup OSp(2n|2n), and the action functional of the WZW model has the 
usual form, 

W[M] = -L [ d^x STr {M-'d^Mf + ^M . (23) 

ioTT Jy, 247r 

Following Witten, the multi- valued functional r[M] is expressed by assuming 
some extension M of M to a 3-ball B that has position space S for its 
boundary {dB = S): 

T[M] = /" STr (M-MM)^3 
Jb 

= f d^x e^^A STr M~^d^M M-^dM M'^dxM . (24) 

Jb 

Of course, in order for this to work we must assume the position space S to 
be a surface without boundary, say a two-sphere or a two-torus. 

The nonabelian bosonization rules tell us to replace the free Dirac the- 
ory plus b-c ghost system by the WZW model with action iy[M], and the 
bilinears 00*^3 and S300* by i~^M resp. i~^M~^. (The factor i~^ is a large 
mass scale which enters for dimensional reasons and depends on the regular- 
ization scheme.) Doing so, and setting T = e^ with X = — S3XS3, so that 
gSg = TJ:3T'^I:s = e^^ = T2, we obtain 

e-^^/^ = fvM exp (-W[M] - ^ f d^xSTr{MT-^ + T^M-')] . 

Next, we remove the factor T^ = e^ from the second term in the exponent, 
by changing integration variables from M to M' = MT'"^: 

s,/N _ f jy^j> g^p f-W[M'T^] -^ I d^x STr(M' + M'-^ 
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(Note that, since the adjoint orbit of T on E3 consists of two disjoint com- 
ponents, this change of variables can only be valid if M runs through two 
components, too.) T^ now appears in the argument of the WZW functional 
W, while the other term has become a plain mass term. 

The last step is to argue that for small length i, or large mass /x, the 
WZW field M' fluctuates only weakly around M' = 1. Indeed, the potential 
STr(M + M~^) has an absolute minimum at M = 1, as follows from the 
definition of the target space as a Riemannian symmetric superspace, see 
equation (p9D below. The leading approximation then is to set M' simply to 
unity, which yields 



'2l 



S,[q] =NW[T 

<7=TS3T-i 

Calculating the first correction due to fluctuations of M' around 1, one finds 
a term with four gradients, multiplied by the square of the length scale 
\Jtj[i. This four-gradient term becomes small, and the approximation of 
setting M' = 1 is therefore justified, on length scales larger than \/ tj \i. 
Note that if the saddle-point approximation is applied after renormalization 
has brought the initially small value of g to about unity, the two length scales 
t and yU~^ ~ 21^ gjix are parametrically the same. It is then this length scale 
[i~^ ~ £ which will set the short-distance cutoff of the nonlinear sigma model. 
By inserting M = T'^ = e^"^ with X = — S3XS3 into the first term of 
iy[M] and comparing with S'eflf[q'] for q = TTi-^T^^ = e^^Es, we readily infer 
/ = 1/A^. This result for / is easy to verify by direct gradient expansion of 
S5, without passing through the WZW model. A more accurate calculation, 
taking into account the finite value of fj,/i, gives 

Now we tackle the delicate task of calculating the coupling 9. Since we 
expect the topological term of the nonlinear sigma model to arise from the 
topological term T[M] of the WZW model, we substitute M = T^ into T[M]. 
Doing so, we naively get zero - the technical reason is that X„ does not 
support such a term - but only naively so. The point to observe is that 
writing r[M] in the form ( P^ requires a smooth extension of M to the ball 
B. If we insist on MS3 = T^Ss = g taking values in X„, such an extension 
does not exist, by the topological obstruction 112 (X„) 7^ (for n > 1). 

We can avoid the obstruction by allowing q to vary over a larger set, say 
the complex group OSp(2n|2n). Then, if < s < 1 is a radial coordinate for 
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the ball B, we can take the extension to be 

M(x, s) = T{x) exp (±is7rS3/2) T{xY^{Ti^^) ■ 

It is seen that for s = 1 we have M(x, 1) = T{x)Ti^T{x)~^T,^ = g(x)E3, 
while for s = we get M{x, 0) = ^iT,^, independent of x. By inserting this 
extension into the expression (|2^) for r[M], and converting the integral over 
B into an integral over E = dB using Stokes' theorem, we find 

i _. „ , 1 
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T[qTj3] = ±— / Sxe^y'^Tiqd^qdyq . 



Given the relation S^[q] = NW[qTj3] and the formula for iy[M], comparison 
of this result with the nonlinear sigma model action (|2l| ) yields 



9 = ±N7r . (26) 

The sign ambiguity arises from the multi-valuedness of r[M]. It has no 
consequence, since the physics of the nonlinear sigma model is periodic in 
the topological angle 9, as long as the position space S has no boundary. 

For a complete description, what remains to be done is to augment the 
effective action with a symmetry-breaking term for finite energy E y^ 0. By 
shifting q \-^ q — i{E /2n)I]s in (p2D, and expanding to linear order in E, we 
obtain the full effective action 

SeIq] = 5eff[g] -'-^Jd'xSTrqE, . (27) 

Finally, a standard calculation starting from the Lagrangian (|^) yields the 
following expression for the local density of states: 

u{E;x) = ^Re [ VqTi qpp{x)a3e-^^^''^ , (28) 

27r5f J 

where the functional integral has to be computed with a positive imaginary 
part Imii^ > 0, in the limit Imii^ -^ 0. 

The above derivation of the couplings / and 9 was only of formal validity, 
as the precise definition of the WZW model was omitted and no justification 
of the bosonization rules was given. In the following two subsections, we are 
going to fill in the gaps. Readers not interested in these details are urged to 
proceed directly to Section p. 
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7.1 WZW model of type C\D 

In the present subsection we will construct the target space of the WZW 
model, based on the general notion of Riemannian symmetric superspace. 
Although this construction is by no means new, it is not as well known 
as it should be. Many influential authors have based their considerations on 
WZW models over the real supergroups GLR(n|n) or OSpi[j(2n|2n), or on the 
unitary supergroups U(n|n). While such practice may have become standard, 
and is admittedly rather convenient and quite satisfactory for a number of 
formal calculations in current algebra, it does not take the WZW model 
seriously as a functional integral, and we are not going to adopt it. Bosonic 
{i.e. non-super) nonlinear sigma models, including the class of bosonic WZW 
models, are defined as functional integrals of maps from a Riemann surface 
into a target space with Riemannian structure. If our purpose is to extend 
this definition to the super world in a mathematically sound way, we have 
to address the fact that the natural or invariant geometry on a supergroup 
invariably is non-Riemann. 

To appreciate the difficulty in some detail, consider U(n|n) for example. 
The even part of the Lie superalgebra of U(n|n) is the direct sum of two 
copies of u(n), which is spanned by the anti-Hermitian nxn matrices. Thus, 
bosonic tangent vectors at the unit element of U(n|n) are pairs {A, B) with 
A^ = —A and B'^ = —B. The natural U(n|n) invariant quadratic form 
(or metric) evaluated on {A, B) is the supertrace Tr A'^A — Tr B'^B, which 
has indefinite sign. As a result, the action functional of the principal chiral 
nonlinear sigma model with target supermanifold U(n|n) is bounded neither 
from below nor from above. On field configurations that fiuctuate rapidly 
in space, the action becomes arbitrarily large, and since the action can be 
either positive or negative, the functional integral is unstable with respect to 
fiuctuations for any choice of sign of the coupling constant. Therefore, unless 
some additional procedure (such as analytic continuation from A^A — B^B 
to A'^A + B'^B) is specified, the theory with target U(n|n) does not exist. 

Such reasoning is not restricted to U(n|n) but holds for the other cases as 
well. Indeed, the Cartan-Killing form of any Lie superalgebra is a supertrace, 
i.e. a difference between two traces. For the supergroups listed above, this 
implies that any rank- two (super) symmetric tensor k, invariant under the left 
and right group actions is necessarily non-Riemann, by which we mean that 
the metric tensor obtained by restriction of k to the bosonic support has 
indefinite signature. For this universal reason, supergroups are ruled out as 
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a target spaces for nonlinear sigma models, at least in the literal sense. 

We understand that the Lagrangians of the "supergroup" WZW mod- 
els which exist in the literature, are used mainly for bookkeeping purposes, 
or as a device to generate equations of motion or other identities, such as 
the Knizhnik-Zamolodchikov equation, that do not depend on working with 
a Riemannian geometry. In the present paper, we wish to put the WZW 
model (or some descendant theory) to a more severe test: we intend to inte- 
grate carefully over the global zero modes (treating the theory in zero-mode 
approximation), so as to establish the connection with exact random-matrix 
limits that have previously been obtained. For that nonperturbative purpose, 
we must construct a functional integral that exists as such. 

While the basic difficulty with supergroup targets is easily recognized, 
it is not so easy to circumvent. For example, switching from U{n\n) to its 
noncompact analog GLc{n\n) /U {n\n) does not improve the situation. The 
latter space is still non-Riemann, and the nonlinear sigma models over it make 



sense only in combination with some procedure of analytic continuation ||26 |. 
It turns out that the difficulty cannot be overcome by using only the tool- 
shed of standard supermanifold theory. What is required is a novel concept, 
namely that of cs-manifold^ which transcends the traditional categories of 
real- analytic and complex- analytic supermanifolds. In short, cs- manifolds 
are super bundles supported by a real-analytic manifold, with a fibre which 
is a Grassmann algebra over C carrying no operation of complex conjuga- 
tion or adjoint. In plain physics language, the "bosons are real while the 
fermions are complex" . It is not hard to see that this is just the right kind 
of mathematical setting to use for the construction of nonlinear sigma mod- 
els with superspace target. Indeed, what we want is a target space with a 
Riemannian metric on its bosonic base manifold, giving an action functional 
bounded from below. A Riemannian metric distinguishes the real numbers 
from the imaginary numbers, in the sense that the tangent vectors with pos- 
itive length form a vector space over R (and not over iR or C). This is what 
is meant by saying that the "bosons are real" . On the other hand, fermionic 
"integration" according to Berezin is nothing but differentiation, and the def- 
inition of Berezin's integral requires no notion of reality, so it is most natural 
to "leave the fermions complex" EH]. Our intention as statistical physicists is 



to do integrals (computing disorder averages of Green functions), and nothing 

"^This terminology is due to Bernstein ||2^. The letter "c" stands for complex, and "s" 
for super. 
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but integrals. Hence, we may take the extreme point of view that complex 
conjugation, while indispensable for the bosons, is a forbidden operation on 
fermions. This is the point of view advocated in [Q, and it is the one we 
adopt here. Note that this means that we bar supergroups such as U(n|n), 
the definition of which requires the use of an adjoint for the fermions. 

Having prepared the stage with these remarks, we now turn to the de- 
scription of the target space of the WZW model. Our starting point is the 
complex supergroup G = OSpc(2n|2n). Its elements, which we denote by 
M, satisfy the equation M~^ = M^, where the orthosymplectic transpose 
M* was defined by (|^) and (M0)* = 0*M*. The group G comes with a rank- 
two supersymmetric tensor k = — STrdM~^dM = STr(M~^dM)^, which 
has the distinctive property of being invariant under left and right transla- 
tions M I— i> gj^Mg~j^^. Let Go = 0(2ra, C) x Sp(2n,C) denote the ordinary 
(or bosonic) subgroup of G. Since Go is a group of matrices with complex 
entries, the geometry induced on Go by restriction of k is non-Riemann. 

What we need to do now is to specify a submanifold of Go on which the 
geometry induced by k is Riemann. This is done as follows. From 0(2n, C) 
(the FF-sector) we select a submanifold Alp isomorphic to the compact or- 
thogonal group 0(2n). Near the group unit, Aip is parametrized by e^ 
where Y = — F^ = —aiY^ai. The group 0(2n) acts on JUp independently 
on the left and right by e^ ^-^ O^e^O]^^, and this action is transitive, which 
is to say that all elements oi Aip are translates of unity. In the BB-sector 
we proceed differently, by selecting the intersection, A^b, of Sp(2n, C) with 
the set of positive Hermitian matrices. The elements of A^b can be written 
as gg"^ with g G Sp(2?7,, C). Since forming the product gg'^ divides out the 
maximal compact subgroup Sp(2n) of Sp{2n, C) on the right, the set A(b 
is isomorphic to Sp{2n, C)/Sp(2n). Using the exponential map, AIb can be 
parametrized by e'''" where X = +X''" = —a2X'^a2- Note that the product 
of two positive Hermitian matrices is not positive Hermitian in general, so 
AIb is not a group. Rather, A^b is a noncompact symmetric space, and the 
complex group Sp(2n, C) acts on AIb transitively by e'^ i— > ge^g^. 

We now look at the product AIb x AIf- This may seem like a some- 
what unnatural hybrid to consider, since Alp is group whereas AIb is not. 
However, both Aip and AIb are Riemannian symmetric spaces, the former 
of compact and the latter of noncompact type, and the product AIb x AIf 
has the desired property of being a Riemannian submanifold of Go. To es- 
tablish the Riemannian property, we first inspect the tangent space at unity. 
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Tl {Mb X Mf). Its elements are pairs A(BB = Aq<a®b) ^ ^^^^ evaluation 
of K = STr {M-^dAiy on these (at M = 1) yields 

STr {A®Bf = TrA'^- Tr B^ = Tr A^A + Tr B^B > . 

From this we clearly see that k, is positive definite on Ti (TWb x -^f)- This 
property carries over to all oi A4b x -^Fi by the invariance of k, under the 
transitive action of the group Sp(2n,C) x (0(2ra)L x 0(2n)R). Thus, k re- 
stricts to a Riemannian structure on A^b x A^f as claimed. 

As an immediate consequence, we have that the numerical part of the 
function STrM = STrM~^ on M.^ x M.p is locally expressed by 

STr Mo = Tr e^ - Tr e^ = Tr e"^ - Tr e'^ 

= Tr cosh y/WI - Tr cos VBW . (29) 

This function has an absolute minimum at M = 1. 

In summary, the object at hand is highly structured, consisting of the 
complex supergroup G = OSp£(2n|2n) with metric k = STr (M~^dM)^, and 
a Riemannian submanifold Ais x A^f- The triple (G, k, Mb x AIf) is what 
is called P] a Riemannian symmetric superspace of type C\D. The name 
encodes the fact that the noncompact BB-sector is symplectic, or C, while 
the compact FF-sector is orthogonal in even dimension, or D. 

The virtue of Riemannian symmetric superspaces, as opposed to super- 
groups, is that they make valid target spaces for nonlinear sigma models, 
as follows. Let the action functional be denoted by S and constructed in 
the usual way, i.e. if K,ij{(f)d(p^dip^ is the expression of the metric in terms 
of target space supercoordinates ip\ which are bosonic for i = l,...,p and 
fermionic for i = p + 1, ...,p + q, we have S = J (Px Kij{ip)df^ip^d^ip^ . Ac- 



cording to Berezin |29[ superintegration is a two-step process. First we do 
the Fermi integral, which is to say we differentiate e~'^ with respect to the 
fermionic fields. (For this no definition of an adjoint or complex conjugation 
of the fermions is needed.) This "integral" always exists - provided we take 
care to control the infinite number of integration variables - in the sense that 
the derivative of an analytic function always exists. The result of doing the 
Fermi integral is a functional, say e~^°DQ, on the bosonic fields. In the sec- 
ond step, we carry out the integral over the bosonic fields, which take values 
in a real target manifold {Mb x A^f in the present case). This functional in- 
tegral exists (modulo the notorious need for regularization in the ultraviolet, 
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regularization of zero modes etc.) because we have constructed the target as 
a Riemannian manifold, with a bosonic action Sq bounded from below. 

One might object that we are violating "supersymmetry" by using com- 
plex conjugation on the bosons (to fix the Riemannian submanifold Mb x 
Aip), while barring the use of complex conjugation and any other adjoint on 
the fermions. This is not so. The role of supersymmetry here is to equip the 
theory with a BRST symmetry, thereby turning it into a kind of topological 
field theory and ensuring normalization of the partition function to unity. 
BRST symmetry does not require that we treat bosons and fermions in an 
egalitarian manner with respect to complex conjugation. The essence of the 
argument can be captured by looking at a simple zero- dimensional example. 
For any analytic function / : M+ U {0} -^ C with rapid decay at infinity and 
/(O) = 1, consider the superintegral 

^ [ dxdyd^d, f{x^ + y' + 2ir,) = /(O) = 1 , 

where x and y (the "bosons") are Cartesian coordinates of R^, and ^ and r] 
(the "fermions") are Grassmann variables. The integral is always equal to 
unity, and this holds true irrespective of whether any conjugation properties 
are imposed on ^ and t] or not. (Thus we do not need to say that S, and r] are 
"real", or rj is the "complex conjugate" of ^, or anything of that sort.) The 
integral reduces to /(O) = 1 by dimensional reduction [^, as a consequence 



of invariance of the integrand under BRST transformations 

(5x = e^ , 6y = er] , S^ = —ey , 5rj = ex . 

Since the generator of the BRST transformation vanishes only at the origin 
X = y = 0, the integrand on M^\ (0, 0) can be written as the BRST derivative 
of another function, causing localization of the integral to the origin. The 
same localization principle applies in the functional setting. 

One might ask again the basic question why we use the highly structured 
concept of Riemannian symmetric superspace instead of more conventional 
constructions (invoking an adjoint for both bosons and fermions). Isn't there 
a simpler way of doing it? We have tried hard, but apparently the answer is 
no, if we insist on global supergroup symmetry and the stringent requirement 
of an invariant Riemannian structure on the bosonic base of the target space. 
There is one thing, however, that we could do in the way of simplification or 
economy of formulation, which would be to eliminate most of the complex 
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supergroup G from the triple {G, k, A^b x -'^f), and retain just the Grass- 
mann algebra fibres over the points of the base M.-Q x M.^- The resulting 
object would be a cs-manifold in the terminology of Bernstein. We do not 
take this step here. 

The WZW model whose existence we postulated earlier for the purpose 
of expanding S^[q\ in gradients, is the functional integral of maps M from 
position space into the Riemannian symmetric superspace of type C\D, 

OSpc(2n|2n),STr(M-MM)2,A^B x Mi 



with the action functional given in (p3D. (To make the model completely 
well-defined, we add an infinitesimal mass term e J(i^xSTr(M + M~^) to 
regularize the zero modes.) We propose to call it the C\D WZW model 
at level k = 1. The simpler name "OSp(2n|2n)" WZW model would be 
inappropriate for two reasons. Firstly, it would belittle the above discussion 
emphasizing the necessity to use a Riemannian symmetric superspace instead 
of a supergroup target. Secondly, such a name could cause some confusion 
since there exists yet another "OSp(2?2|2n)" WZW model, still at level fc = 1, 
which differs from the present one by the exchange of the BB- and FF-sectors 
- in other words, the bosons carry an orthogonal structure while the fermions 
are symplectic. (The latter model, with k = 2, may have some relation to 



the fixed point governing the spin quantum Hall transition |^ ^ Q.) 

There is one point that needs further attention. Although the existence 
and stability of the principal chiral nonlinear sigma model with target C\D 
should now be clear, the situation for the WZW model is still precarious. 
The reason is the presence of the multi- valued term r[M] in addition to the 
kinetic term induced by the metric tensor n. While ir[M]/247r is imaginary 
(and ambiguous by integer multiples of 27ri) in the FF-sector, it is real (and 
single- valued) in the BB-sector. Indeed, ii A = A^ , B = B"^ and C = C^ are 
elements of the tangent space of A^b at unity, the 3-linear form iTT{A[B, C]), 
which induces the BB-part of zr[M], takes values in M. Moreover, this term 
reverses its sign under {A,B,C) \-^ {—A,—B,—C), so there exist field di- 
rections along which q-^^I^]/'^^^ increases exponentially, thereby jeopardizing 
the existence of the functional integral with action iy[M]. Fortunately, one 
can prove the following bound: 



Re 

247r ^ ^ 



<^jd'xSTi{M-%M)\ (30) 
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which ensures the stabihty of the functional integral Je"'^'*^'. For the 
special case A^b = H'^ (three- hyperboloid) , which is isomorphic to TWb = 
Sp(2n, C)/Sp(2n) for n = 1, this bound was established in p^ . 

As a final remark, let us mention that the bound ( pO]) is optimal. In other 
words, if we multiply the kinetic term of the WZW action IV [M] by a factor 
less than unity, and arbitrarily close to unity, there exist field configurations 
violating the bound, and the functional integral becomes unstable. This 
means that the C\D WZW model resides exactly on the border of stability. 
An important consequence is that iy[M] tolerates the addition of a JlJr 
current-current perturbation only with a definite sign of the coupling. 

7.2 Nonabelian bosonization 

With a good functional integral in hand, we can now go ahead and establish a 
supersymmetric extension of nonabelian bosonization. We are going to show 
that the free Dirac theory plus b-c ghost system is equivalent to the C\D 
WZW model at level k = 1. To that end, we consider the partition function 

Z^zwlA A]= f VM exp -II^[M; A, A] , 

where IV [M; A, A] is a gauged WZW action: 

W[M;A,A] = W[M]-— I d^xSTrfM-^dMA + AMdM-^ 

+Amam-^ - AA^ 

and A and A are external sources taking values in osp(2n|2n). Our main tool 
for analysing the WZW model is the Polyakov-Wiegmann relation |35 

W[gh~^] = W[g] + W[h-^] - — [ d^x STr g-^dg h^^dh . (31) 

2n J 

By using it, we easily derive the identity 

W[g,Mg],'] = W[M;gl'dg^,g],'dg^] + W[g,g^'] , 
and parametrizing the external sources by 

A = gl'dgL, A = g],'dg^, (32) 
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we can compute the partition function Zwzw[^)^] exactly: 

Zwzw[^,^] = I VM e-^^3LMg-a'\+w[g^a-^'] = expiy[^^^^^] . 

To arrive at the last expression, we changed variables from M to M' = 
g^Mg^ in the functional integral. Such a substitution is certainly valid at 
the infinitesimal level, i.e. as long as the Taylor expansion of ^vvzw[^5^] 
with respect to the sources A, A is truncated at finite order. 

From the Polyakov-Wiegmann relation, one can also verify the invariance 
of PF[M] under local 0'&\i^{2n\2n) l x OSpc(2n|2n)ij transformations. 



M{z,z)^gL{z)M{z,z)gR{z 



^\-i 
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which is characteristic of any WZW model and takes the form of an afiine 
Lie symmetry in the quantum theory, which here is nonunitary. (Again, 
the invariance holds without doubt at the infinitesimal level, since we can 
always pass to the complexified tangent space with impunity.) By Noether's 
theorem, the invariance entails two sets of conserved currents: 

J = MdM-^ , J = M-^BM , 

satisfying the equations of motion dJ = dJ = 0. 

Consider now the free Dirac theory plus b-c ghost system, with the par- 
tition function 

This functional (super) determinant is ill-defined and needs to be regularized 
in the infrared, say by adding to the Dirac Lagrangian a small mass term, 

e (0*S30 + 0*S30) = 2e(|6p + \c\^ + ...). 

Using standard heat kernel techniques [^, the regularized superdeterminant 
can then be computed to be the exponential of another WZW action: 

ZDirac[A, A] = SDet-1/2 (^ ^'^'^ ^ +^^ ] = exp Wlg^Jlsgn'^s] , 
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where the sources A and A were again assumed to be of the form (3^). By 
comparing answers, we see that the partition functions of the C\D WZWfc=i 
model (derived from the complex supergroup OSp(2?T,|2?T,)) and the Dirac+ 
ghost system (with n species of particles) coincide: 



on gauge source fields of the form (^2]) . This equivalence furnishes the basis 
for nonabelian bosonization. 

As an immediate consequence we get bosonization rules for the currents: 

00t ^ —dMM'^ , S300*S3 -M-^BM , (33) 

2n 2% 

by comparing the J ■ A perturbations in the two theories. (Using these 
bosonization rules we can rewrite a single species of Dirac fermions with 
random mass as a WZW model with marginally irrelevant current-current 
interaction. No additional insight seems to be gained from that alternative 
representation.) However, in our computation of the gradient expansion of 
5*5 [g], Eq. (^21), we employed the bosonization rules 



00*^3 ~ {27ii)-^M , and ^300' ~ {2ni)-^M-^ . (34) 

To justify these, further considerations are necessary. 

The C\D WZWfc=i model, as any WZW model, is solvable to a high 
extent, owing to the availability of exact operator product expansions be- 
tween the conserved currents and the fundamental field M (and all other 
primary fields). To write these down, we introduce the convenient nota- 
tion J^ = —^STiAMdM'^, where A is any spatially constant element 
A G osp(2|2). Then, making in the functional integral {M{w,w)) a change 
of variables M -^ e~^^M, where X is holomorphic in a small neighborhood 
of the point w and smoothly goes to zero outside, we obtain the operator 
product expansion |5^ 



J'^(z)M(w, w) = M(w, w) + ... . (35) 

z — w 

For the antiholomorphic current defined by J^{z) = —^STiAM'^dM, an 
analogous calculation gives 

J^(z)M(w, w) = M(w, w)- + ... . (36) 

z — w 
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The OPE of the conserved currents among themselves can be obtained in a 
similar manner. For the holomorphic current we have 

,.(.„>,) ^_Mli^,Z^^,..., (37) 

[z — w)^ z — w 

and the same formula (with z —>■ z etc.) holds for the antiholomorphic one. 
On general grounds, the holomorphic component T of the stress-energy 
tensor has the canonical Sugawara form T{z) oc Hij : J^{z)J^{z) : , where J* = 
— |STre*M(?M~^, and {e*} is a basis of the Lie superalgebra osp{2n,2n). 
The metric tensor is expressed by Hij = — STrejCj, indices being lowered 
via 6j = STre*ej. Classical considerations (by the construction of T from 
the Lagrangian via Legendre transform) would suggest a constant of propor- 
tionality 1/k, where /c = 1 in the present case. However, the constant is 
renormalized by quantum fluctuations. Its correct value is deduced from the 
requirement 



tA 



JHw] 



Tiz)J^iw) = -^^ + ..., 
[z — wY 

expressing the fact that J is a holomorphic conserved current and therefore 
must have conformal dimensions (A, A) = (1,0). Expanding the product 
T{z)J^{w) with the help of (^) and the associativity of the operator product 
algebra, one finds 

T{z) = -^ : J\z)J^{z) : , (38) 

where h^,, called the dual Coxeter number, is the quadratic Casimir invariant 
evaluated in the adjoint representation: 

Kij[e\[e^,A]] = KA. 

The superbracket [■, ■] here means the commutator or the anticommutator, as 
is appropriate. A standard calculation using the root system of osp(2n, 2n) 
yields 

/i* = -2 , 

independent of n. Note that /i* is negative! We will see shortly that this has 
to be so, in order for the dimensions of operators such as the fundamental 
field M to come out being positive. (We can also check the sign of /i* by 
interpreting the C\D WZW model with say, n = 1, as the replica limit r — > 
of the familiar 0(r) WZW model. As is well-known, the Lie algebra of 0(r) 
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has dual Coxeter number h^, = r — 2, which indeed becomes —2 at r = 0.) 
The Fourier modes of T{z), 



(27ri)~i (p T{z) z'^^^dz 



are the generators of a Virasoro algebra, with the central charge c being zero 
by supersymmetry. Similar statements hold for the antiholomorphic compo- 
nent T{z) of the stress-energy tensor. The presence of these Virasoro algebras 



implies the conformal invariance of the WZW model. From equations ([351)- 
( PSP one sees that the leading singularities in the operator product expansion 
of the stress-energy tensor with the fundamental field M are 

T{z)M{w,w) = . ^^\ m{w,w) + ..., 

[z — wy 



[z — w)2 



where 



A - A ^^ 



and Cm = HijC^e^ is the quadratic Casimir in the fundamental representation 
of osp(2n|2n). The numbers (A^, ^m) are the conformal dimensions of the 
field M. We can calculate Cm by making an explicit choice of basis {e*}. In 
this way we find Cm = ~l/2, which yields a total dimension of one: 

dim(M) = Am + Am = 2 x -^^ = 2 x ^^ = 1 . 

This means that the two-point function (MM) must vary as the inverse 
square of the distance. Global Gl x Gpt invariance then constrains the OPE 
of M with itself to be of the form 

STi AM(z,z) X STi BM-\w,w) = -- — ■ —STtAB + ... , 

{2n\z — w\)^ 

where £ is a length scale depending on the choice of UV regularization scheme. 
The supermatrices A and B are arbitrary elements of osp(2n|2n). 

We now return to the Dirac-|-ghost system, and rewrite the OPE given 



m 



11) in the invariant form 



STr AB 

STr A0(^)0*(^)S3 X STt B^,^{w)<l)\w) = --— ^ + ... 
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(The presence of the factors S3 leads to the overall minus sign on the right- 
hand side.) By comparing the two OPEs we see that they agree to leading 
order if we identify M/ d, with 00*S3, and M~^ ji with S30(/)*. Moreover, the 
OPEs between these fields and the respective currents agree to leading order 
in the two theories. By computing the three- and four-point functions one 
finds [pTl that the agreement persists to higher order. This then justifies the 



bosonization rules (|3^). Their status is the same as in the usual case |^ 



In summary, we have established the C\D WZWa:=i model to be equiva- 
lent to the Dirac plus h-c ghost system at the current algebra (or infinitesimal) 
level. What about aspects that transcend the infinitesimal level? This is not 
an empty question, as the Riemannian symmetric superspace of type C\D 
with bosonic submanifold M^ x M.^ = (Sp(2n, C)/Sp(2n)) x 0(2n) has 
two connected components, owing to the topological equivalence 0{2n) ~ 
Z2 X S0(2n). Since the discreteness of Z2 allows the WZW field to break up 
into domains separated by domain walls, we are led to ask: exactly how is 
the functional integral to be defined? Do we have to sum over all possible 
numbers and positions of the domain walls, or do we not? 

We believe that the answer to the question is yes, and the precise non- 
perturbative definition of the functional integral does have to involve Z2 as 
a local degree of freedom, for the following reason. Z2 already exists in the 
classical setting, and is not a peculiarity of the supersymmetric formalism. 
Indeed, according to Witten |jl^, the free 2d massless n-component Majorana 



fermion theory is equivalent to the level-one WZW model not over SO(n) but 
over the disconnected group 0{n). The argument given in [Q proceeds by 
the comparison of current algebras, and does not address the issue of do- 
mains and domain walls in the Z2 degree of freedom of 0{n). However, by 
specializing to the case n = 1, we see that Z2 does need to be summed over 
locally, for a single Majorana fermion is not an empty theory. Rather, its 
partition function bosonizes to that of the 2d Ising model, which is a theory 
of local Z2 spin degrees of freedom. In the continuum limit near critical- 
ity, the latter partition function is computed by summing over all possible 
numbers of domain walls and their positions. By counting the number of 
degrees of freedom we see that the situation is the same for n > 1: the 
fermionic Hilbert space contains states that cannot be produced by acting 
with currents on the vacuum. To achieve equality of the partition functions, 
one needs Z2 in the WZW model, and Z2 must be summed over as a local 
spin degree of freedom. We particularly emphasize this point, although we 
will not pursue it seriously in the present paper. 
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Here ends our extensive excursion into nonabelian bosonization and the 
C\D WZW model. This material was included to make the present paper 
self-contained. 

8 Density of states 

In Section ^ we saw that, when fluctuations are neglected, the composite 
field Q assumes a nonzero saddle-point value /x. Since Q enters the theory by 
coupling to the Dirac bilinear (j)i(f)\ + 0^0;, the latter acquires a nonvanishing 
expectation value, too, and the density of states at zero energy becomes 
finite. This happens in spite of the fact that the Dirac theory is devoid of 
any scale, and is an instance of dynamical mass generation. Of course, a 
fixed value of the "order parameter" Q breaks the OSp(2rz|2n) symmetry 
and leads to the existence of Goldstone modes, the low-energy effective field 
theory for which is the nonlinear sigma model (pi]). 

What is the effect of order parameter fluctuations? The Mermin-Wagner- 
Coleman theorem states that continuous symmetries cannot be spontaneous- 
ly broken in two dimensions. This statement applies to compact symmetries, 
or nonlinear sigma models with compact target spaces. These flow under 
renormalization to strong coupling, so that the field fluctuations grow large, 
and the potentially broken symmetry is restored in the infrared. However, 
for nonlinear sigma models with supersymmetry, or zero replica number, 
another scenario is possible, as has long been known from the example of 
time-reversal invariant disordered 2d electron systems with spin-orbit scat- 
tering (class All). In that case, the beta function at weak coupling has the 
"wrong" sign, which physically corresponds to weak anMocalization, and 
the nonlinear sigma model undergoes logarithmic flow to a Gaussian fixed 



point describing a perfect metal. As was already mentioned in [^ , the same 



happens in the present case. Let us review how to arrive at this result. 

8.1 Renormalization group 

Friedan HO] has shown in great generality that the one-loop beta function 



of a nonlinear model with target space metric k is determined by the Ricci 

curvature R: 

dn R 

dlni ~ ~2^^ '" ■ 
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The metric of the present target space is k = — (87r/)^^STr (dg)^, and its 
Ricci curvature can be described as follows. (Note that the topological term 
of the nonlinear sigma model does not renormalize at weak coupling, and 
can safely be ignored. We shall also neglect the nonperturbative Z2 degree 
of freedom of the target space.) We parametrize q by the exponential map, 
q = e^Ti^e'^ with S3XS3 = —X, and denote the commutator (or adjoint) 
action by [X, •] = ad(X). Then the Ricci curvature is the second-rank 



invariant tensor determined by the quadratic form [^ 

Ro{X,X') = -STrad(X)ad(X') 

on the tangent space at X = 0, or g = S3. For any irreducible symmetric 
space, this tensor will always be a constant multiple of the metric tensor, 
as follows on general differential geometric grounds and is necessary for the 
nonlinear sigma model to be renormalizable. Moreover, the constant of pro- 
portionality is independent of n by supersymmetry, and we may calculate it 
by looking at the simplest case n = 1. To do so, one evaluates both tensors 
on some element of the tangent space, say H = HE-qb ® en- The eigenvalues 
of the adjoint action are called roots. In the present case with n = 1, there 
exist two nonvanishing roots, a bosonic root 2h with multiplicity one, and a 
fermionic root h with multiplicity two. Hence, 

-Ro{H,H) = {2hf - 2h^ = 2h^ = STtH^ , 



which gives R = STr (dg)^/4. In terms of the coupling /, the one-loop RG 
equation then reads 

-f . (39) 



dlni 
By integrating this equation from the cutoff scale £0 to i, with initial condition 

/o = /(£o) = X-i(l-e--/^)-\ 

we find that the coupling / flows to zero as an inverse logarithm: 

This results in the system being a perfect "metal" , with a resistivity (~ /) 
that goes to zero in the thermodynamic limit. By the same token, the broken 
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OSp(2n|2?2) symmetry will not be restored, but remains truly broken, and 
the density of states at i? = is sure to be nonzero. 

As was mentioned earlier, spontaneous breaking of continous OSp(2n|2n) 
symmetry also occurs for 2d metals with spin-orbit scattering (class All). 
There is, however, one important difference we wish to mention. For class 
y4II a metallic phase does not exist in disordered wires (the quasi-ld limit), 
where weak antilocalization at short scales always crosses over to strong 
localization at large scales. In contrast, for class D anomalous "metallic" be- 
havior already occurs in quasi-ld systems! By solving the CI | Dili nonlinear 
sigma model using its quantum Hamiltonian, one finds that the (thermal) 
conductance decays algebraically for wires of arbitrary length. The anoma- 
lous behavior can also be seen on a more elementary level from the maximum 
entropy transfer matrix ensemble, evolving according to the so-called DMPK 
equation, for class D.^ 

8.2 RG for the density of states 



We now embark on a quantitative calculation, based on (|28D , of the local 
density of states in the metallic limit. The idea is to follow the flow of 
the renormalization group starting from the short- distance cutoff io all the 



way up to the system size L, and then evaluate the functional integral (2S) 



in zero-dimensional approximation. (The local Z2 degree of freedom of the 
nonlinear sigma model is expected to play no essential role in that process, 
and will neglected but for its global part.) The initial cutoff is taken to 
be the length associated with the dynamically generated mass, io ~ fi~^, 
which sets the scale over which the "ballistic" free fermion theory crosses 
over to the nonlinear sigma model describing diffusion. (From the nonabelian 
bosonization argument given at the beginning of Section ^ we saw that the 
reduction to the nonlinear sigma model is valid on length scales larger than 

To begin, we add a source term A Tr qpp{0)a3 to the action functional and 
differentiate with respect to A at A = 0. From [^, the parameter A obeys 
the renormalization group equation 



din 



^We thank John Chalker for pointing this out. See also 
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where Cq is the quadratic Casiniir invariant of osp(2n|2n), normahzed ac- 
cording to the metric — STr (dg)^/8 and evaluated on the representation the 
field q transforms under. This number turns out to be negative: 

so that A grows with increasing scale. In the thermodynamic limit, this will 
lead to a logarithmically divergent density of states, as follows. 

Let the system be finite and of size L. By integrating the fiow equation 
for A from the microscopic cutoff Cq up to L, 

A(L) = A(£o)exp / /(f)rflnf , 



and inserting the scale dependence ( ^OD of the running coupling /, we obtain 

A(L)/A(£o) = l + /oln(L/4). 

To keep track of the renormalization of the operator j d?x STr gSa in the 
action functional Se-, it is convenient to introduce a running energy parameter 
e(£), with initial value eo = fiNE/2g. This variable, as compared to A, 
carries two extra dimensions from d'^x, and therefore evolves according to 
the equation 

which integrates to 

e(L)/eo = (L/£o)'(l + /oln(L/^o)) • 

Given (|28|) , all this leads to a scaling relation: 

uiE, fo)L = (1 + /o HL/io)) y {Et{L)/to, f{L)),^ . (41) 

Here i'{E, /o)l is the local density of states of a system of size L, energy 
E, and nonlinear sigma model coupling /q, with the short-distance cutoff 
being understood to be Cq. For L ^ Cq and i^^ = 0, we immediately obtain a 
logarithmic law, 

KO,/o)L~/oln(L/£o), 



which is a result given before by Senthil and Fisher |TT 
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Actually, we can work out a more precise answer. The function z/(.)^g on 
the right-hand side of (^l]) is meant to be evaluated for a system whose size 
L equals the cutoff £o- Under such circumstances, the functional integral can 
be calculated by retaining only the spatially homogeneous mode q{x) = go 
(zero- mode approximation): 

i^iEh, = ^ Re I Dqo Tr{q^^as) exp (^^^^ STr q^J^s) , (42) 

which no longer depends on the coupling /. The computation of this integral 
is the subject of the next section. Using the answer given below in (|46|), we 
have 

.(E,hV^H+ \^^^, > . (43) 

where 

p^ = i^{l + f,\n{L/eo)) . 
■ng 

Since we have neglected the influence of finite E on the RG flow, the validity 
of this result is restricted to small energies E <^ Et^\^. By standard reasoning, 
the relevant energy scale is the Thouless energy i^xh = D / L?', where D = 
g/{2niJ,Nfo) has the meaning of a diffusion constant. In the opposite regime 
E S> -E^Th, direct use of one-loop perturbation theory yields 

i^ r cPk 1 



v{E) = Ui^ + -Re / 



(27r)2 Dk^ - 2iE 

where the momentum integral was cut off by |A;| < l/io- We observe that 
the energy- dependent part of the density of states behaves as \ia{l/E) on 
intermediate scales, and as E'"^ in the asymptotic regime of large E. 

8.3 Random matrix limit 

For the case n = 1, we are now going to compute the superintegral 



Ke) 



i / DgTr(gppE3)expiieSTrgE3 , 
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where Dq is the invariant Berezin measure on the Riemannian symmetric 
superspace Xi (Section |]) normahzed by /-Dg exp ^ieSTrgEa = 1. This in- 
tegral appeared in (^), and can be viewed as the density of states in the hmit 
of small system size (also referred to as the ergodic regime, or the univer- 
sal random matrix limit). Although J(e) can be shown to be semiclassically 
exact, which is to say that /(e) is calculated exactly in saddle-point approx- 
imation (by a supersymmetric generalization of the Duistermaat-Heckman 
theorem [0), we shall not make any use of this deep fact and proceed in a 
more pedestrian way. 

The first step is to write down an explicit parametrization for the super- 
matrix q. To prepare that step, it is convenient to change the arrangement 
of the multiplets 0, 0* to 



c 
V b J 



0* = (^+,6,V^_,-c) . 



Next recall the definition of the orthosymplectic transpose of a supermatrix 
X by (X0)* = 0*X*. Elements of the Lie algebra osp(2|2) are solutions of 
the equation X = —X*, and are easily verified to be of the form 



X 



(a 

6 


V7 



a 
d 

7 
c 





P 
—a 

a 



b 

-S 

-dj 



(45) 



where Roman letters stand for commuting numbers, while Greek letters de- 
note Grassmann variables. 

Now, as was emphasized in Section H, the superspace Xi consists of two 
disjoint pieces. The first contains the diagonal matrix S3 G osp(2|2), which 
rearranges to 

S3 = diag(l, 1,-1,-1). 

On this piece of Xi, the supermatrix q will be parametrized as 



Q 



(1) 



S3 + 2 



f ZZ{l-ZZ) 

I z{i-zz)- 



-1 



-Z(l- 
-ZZil 



zz)-^ 
- zz)-^ 



This can be seen to obey the nonlinear constraint g^ = 1. From ( ^5] ) and the 
fact that q = TTj^T~^ lies in osp(2|2), the 2x2 supermatrices Z and Z have 
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to be of the form 



-/? 

13 b 



7 
7 c 



Following Section P, the bosonic base of Xi is fixed by the conditions 

b = c*, |6p<l. 

The second piece of Xi contains the diagonal matrix diag(— 1, 1, 1, — 1), 
which can be represented as 



diag(-l, 1,1,-1) = 0^30-1 



with 



= EFF(S)cri + Ebb ® h 



/O 1 0\ 

10 

10 

\0 1/ 



G OSp(2|2) . 



Since O lies in the orthosymplectic Lie supergroup, the q matrix for the 
second piece of Xi can be obtained from the first one by conjugation with 
O: 

q{2) ^ T'O^^Q-^T''^ = 0Ti:3T-^0-^ = Oq^^^Q-^ . 

The next step is express the integration measure Dq in terms of Z and 
Z. One way of doing it is to use the coordinate expression of the metric 
STr (dg)^. A straightforward calculation yields 

-STr {dqY/8 = STr (1 - ZZ)-^dZ{l - ZZy^dZ . 

This now gives rise to a Berezin measure in the standard way, cf. Appendix 



F of Ref. [0], where a similar calculation was described in full detail for 
symmetry class C. One finds 



Dq = D{Z, Z) SDet (1 - ZZ) 



-1 



where D{Z^ Z) is a fiat Berezin measure. On setting h = re**^, we arrive at 
the coordinate expression 

Dq = (27r)-Vdr A dipdfsd^ o (1 - r^ - 2/^7)-^ , 
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which has been normahzed so as to satisfy / Dq exp |ze STr gSa = 1. 
We are finally ready to compute /(e). Using 

jSTrg^^Sa = - + 



4 ^ ' l-r2 ' (l-r2)2 

1 2r2/37 



iSTrg(2)s 



3 



4 ^ ■^ 1 _ ^2 ' n _ ^2A2 ' 

iTrgJp^j^cJa = 1 - z n = -^Tr q^pas 



and making the substitution t = (1— r^) ^, we obtain the following expression 
for the integral: 



°° rft . . ,. ^2 



/(e) = J ^dpd,{t + 2t'f3^){l-2tf3^)x 

X ('exp (ie (t - 1 + 21"^ (3-/)) - exp (k (t + 2t(t - l)/?7))) . 

where we have combined the contributions from the two pieces of Xi. The 
integral is now easily calculated to be 

1 e" 
/W^l-- + -. (46) 

On taking the real part, we get 

^ ^/ N sine 
Re/ e =1 + , 



in agreement with the known result for the density of states (in scaled 
units) of the Gaussian random matrix ensemble for class D. 

Let us remark that exactly the same result for /(e) would have been 
obtained by treating the superintegral in saddle-point approximation, in the 
spirit of the Duistermaat-Heckman theorem. The constant part comes from 
a saddle point on the trivial component (with respect to Z2) of the target 
space, and the oscillatory part e*'^/ze from a saddle point on the nontrivial 
component. 
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9 Disordered d-wave superconductor 

An important physical realization of symmetry class D is by the low-energy 
quasiparticles of dirty superconductors with broken spin-rotation and time- 
reversal symmetry. Although there is a number of interesting cases to con- 
sider [0, |l3l, we will here focus on a specific model of a disordered d-wave 
superconductor. The treatment in this section follows Ref. [|5] to some ex- 
tent. We shall work in two-dimensional space with Cartesian coordinates x, y 
and wave vector k = {k^, ky). 

The effective quasiparticle Hamiltonian of the pure system is 

^0 = Y^^Sk - fi) clc,^ + Y, Mc^cl,^ + h.c.) , (47) 



where the single-particle energies Ek, shifted by the chemical potential /i, and 
the gap function A^ are taken to be 



= —t (cos k^a + cos kyo) , 
= — Aq (cos k^a — cos kya) , 



Afc = 

and a is the lattice constant of the tight-binding model that results on trans- 
forming to position space. The excitations of this Hamiltonian are gapped 
almost everywhere in the wave vector plane, the exceptional places being the 
four points given by \kx\ = \ky\ = a~^ arccos(— /i/2t), where e^ — n = A^ = 0. 
Since the low-temperature behavior of the system on large distance scales is 
expected to be dominated by the quasiparticles with low energy, it is natural 
to linearize Hq around the nodal points kr^^y = ±a~^ arccos(— yu/2t), denoted 
in a self-explanatory notation by (±±). By following a standard procedure 



18, ^, pOf, we obtain four Majorana Hamiltonians, one for each node. The 



one for (++) reads 

Ht^^ = \^Vl I d'x J2 {fldJl + fadlfa) 

with the two velocities given by vi = at and V2 = ctAo, and we have intro- 
duced rotated coordinates by xi = |(a; + y) and X2 = ^{x — y). The fermion 
operators /„ and fl are linear combinations of the operators c^ and cj. in 
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The Hamiltonian Hq for the nodal point k^ = —ky = 7r/2a is obtained 
from Hq by exchanging di <-> d2- The remaining ones, Hq and Hq , 
are gotten from the previous two by reversing the overall sign. 

We now introduce disorder into the problem. In the present context, 
disorder can be classified into two categories: mtranode scattering, and in- 
temode scattering. We will deal with these scattering mechanisms separately. 
First we take into account the intranode scattering, and later we will include 
the internode scattering as a perturbation. Such a two-step procedure makes 
sense if the "backscattering" between nodes is weak compared to the "for- 
ward" scattering within a node. 

In the family of symmetry classes of disordered superconductors, class 
D is distinguished by the absence of any special symmetries such as time- 
reversal or spin-rotation invariance. Thus we are to add disorder of generic 
type. Without much loss, we can restrict the disorder to be local (involving 
no derivatives). To write down the most general, local, intranode scatter- 
ing Hamiltonian, we again single out the node (++) as an example. The 
expression for this Hamiltonian in terms of the operators /^, fl reads 

where A now is a complex "random gap function" , and Mab = M^a is a 2 x 2 
"mass matrix" fluctuating randomly in space. They are given as certain 
linear combinations of the perturbations (random scalar potential, random 
order parameter, etc.) that are added to the original pure Hamiltonian Hq. 

Consider now the problem posed by the sum H^'^'^^ = Hq + iJJ 
This is essentially the problem we solved in Sections ^^. To make the 
correspondence precise, we choose anisotropic units of length in the xi and 
X2 directions such that f i = f 2 = 1 (we will restore the proper length units 
shortly), and write 

H^"-"-^ = ^/rf^x(/i/J-7^(++)-(^j)+const, 



( Mn M12 %d^^d2 A \ 

M21 M22 -A z^i + ^s 

idx - <92 -A -Mil -M21 

\ A %dx-d2 -M12 -M22 / 



All information about the second-quantized Hamiltonian H^'^^^ is encoded 
in the first-quantized Hamiltonian 7i*^"''+\ and we can equivalently work with 
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the latter instead of the former. Doing so, and adopting white-noise disorder 

(M„b(x)M,rf(x')> = g 5ac5hd 5(x - X ) , 
(A(x)A(x')> = ^5(x-x'), 

with (Mab) = (A) = 0, we arrive at the N = 2 version of the random 
Hamiltonian treated earher. Drawing on the results of Section 0, we can 
immediately write down the low-energy effective field theory for this problem. 
(Of course, the anisotropic choice of length units has made the short- distance 
cutoff of the field theory anisotropic. This does not make a big difference, 
as the nonlinear sigma model couplings are essentially cutoff independent.) 
On putting the proper units of length back in place, we get an anisotropic 
theory (g = g'-"*"'"''), 

5^++) = — [ d^xSTT(-diqdiq+-d2qd2q] 
Svr J \V2 Vi J 

+ — d^xef,^STiqdf,qd^q. 

Note that the topological term by its nature is ignorant of all length units and 
therefore cannot depend on the ratio f 1/^2. Exactly the same effective action 

governs the field q^ •* for the node ( ). The actions for the remaining 

two nodes are obtained by interchanging di ^->- 82, which has the particular 
consequence of reversing the sign of the topological coupling. 

Finally, we turn on the scattering between nodes, which couples the four 
theories together. Here we can easily guess without any calculation what is 
going to happen. The fields q^^^^ for the nodes (st) are the Goldstone modes of 
a broken orthosymplectic symmetry. In the absence of internode scattering, 
the four sectors described by the g*^**-* are decoupled, and the symmetry group 
consists of four independent copies of OSp(2n|2n). Scattering between the 
nodes reduces the symmetry to the diagonal subgroup, which acts by the 
same factor in each copy. On these grounds, we expect that the effect of 
internode scattering is to "lock" the fields of the four nodes to each other: 
q^'^'^^ = q^ ) = q^~^~^ = q^~^^ = q, at large distance scales. The locked field 
q is the Goldstone field due to the remaining diagonal OSp(2n|2n) symmetry. 
The effective action for it will simply be a sum of actions: 

cfx STr dnq dnq . 
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In this final expression the anisotropy of the single- node actions has cancelled, 
and ^eff is isotropic. Also, the four topological terms have added up to zero. 
This was to be expected, since the presence of these terms violates parity, 
which is a good symmetry of the superconductor, unless an orientation is 
induced by a strong magnetic field or by an order parameter with nonzero 
chirality (such as dx2_y2 + id^y)- 

Thus we have arrived at the nonlinear sigma model for class D. The 
model is at weak coupling if the ratio fi/f2 is large, which is the typical 
situation in experimental systems. The coupling f~^ = A^vf + f|)/fif2 or 
rather, const x Tf~^ where T is the temperature, has an interpretation as 
the thermal conductivity of the disordered superconductor [ll|. (A nonzero 



topological angle 9 would have corresponded to a thermal Hall conductivity.) 
The renormalization group for the weakly coupled theory was worked 
out in Section |^. We found that the theory flows to a Gaussian fixed point 
describing a perfect (thermal) metal, and the local density of states in the 
thermodynamic limit diverges logarithmically at ii^ = 0. 

10 Free fermions and phase diagram 

We have shown that a model of N species of Dirac fermions in class D sup- 
ports a metallic phase in two dimensions. Our analysis has put the existence 
of that phase on solid ground for iV ^ 1, or A^ = 2 with large and oppo- 
site anisotropics in the velocities. Here and in the next section, we wish to 
go further and address two questions that emerge from the recent literature 
P, |Tl|, |13[: (i) Can the 2d metallic phase of class D also be realized in the 



constrained parameter space of the fundamental case A^ = 1? (ii) What is 
the location of the free-fermion point relative to the metallic phase? 

The answer to the second question is that the metallic phase has the free- 
fermion point sitting right on its boundary. We are going to demonstrate 
that this is so, by employing the supersymmetric extension of nonabelian 
bosonization developed in Section 0. As for the first question, we will see 
that the answer is no, but there exists a remarkable twist to the story. 

We start with the second question. Recall the Lagrangian L2 in (p!0D, 
describing A^ species of massless Dirac fermions plus b-c ghost system, weakly 
perturbed by the operators $*-°'' with couplings ga (« = 1, •••,4). As we have 
seen, a single species in the pure limit is equivalent to the C\D WZWfc=i 
model with field M and action functional 14^ [M]. We now bosonize L2, by 
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exploiting the equivalence for each species separately: 

ct>\d<Pi + ^\d^,^W[Mi\ (/ = l,...,iV). 

This scheme can be justified by introducing auxiliary fields Q as in (JT^) to 
transform L2 into a Lagrangian for A^ decoupled species. The perturbations 
$*^°) are then bosonized using the rules derived in Subsection |7.2| , which 
results in 

N / N 

S'. = E(^[^']-^7^§^/^'-STrJ,S3J.S3j 
^ E / ^'^ (I STr M.EsM.Es + | STr M^^M, 

+ |(STrMfcE3)(STrM;S3)-^STrJ,S3J/S3) . 



+ 



N 



Here we isolated the part of the perturbation that acts within a single species 
/, and bosonized it by reduction to the bosonization rule for the currents: 

{^\ct>,f = ~STT{Ml)iM\) ^ -(27r)-2STr J,S3J,S3 . 

We observe that, since the fields J; = M^OMj^^ , Ji = Mf^dMi and M; have 
conformal dimensions (1,0), (0, 1), and (1/2, 1/2), respectively, bosonization 
has preserved the marginality of the perturbation. Of course, the couplings 
will not be truly marginal but will evolve under renormalization in a way 
that is determined by the operator product expansion among the various 
perturbations. Given that the perturbed C\D WZWfc=i model is a true 
image of the original theory, the fiow equations are identical to the ones 



worked out earlier and given in ([121) . Note also that, by taking the beta 
functions from the Dirac representation, instead of recalculating them in the 
bosonized theory, we avoid the nontrivial question of how to renormalize the 
local Z2 degree of freedom of the perturbed WZW model. 

To go further, we distinguish between cases. Let us first assume N > 2 
and take the bare couplings to be 5^1(^0) = 5*2 (^o), and 5^3(^0) = 94(^0) = 0. 
Numerical integration of the one-loop fiow equations ([12]) then shows that 
the couplings gi,g2 grow {gi more strongly than (72), fi's becomes nonzero and 
positive, and g^ moves to negative values but respects the bound 5^1 + 5'4 > 
dictated by Hermiticity of the random Hamiltonian. (For N = 2, the fiow 
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eventually is attracted to the line Qi + Qa = 92 = 93 = 0, which has a higher 
symmetry, namely that of class All; see Section ^. Thus the symmetry is 
dynamically enhanced in that case and the flow, though starting from a 
point in class D, terminates in class All.) Hence we are facing a relevant 
perturbation that drives the system away from the free-fermion point. We 
make the reasonable assumption that the relevant nature of the flow persists 
beyond the one-loop approximation used in deriving the flow equations. 

The fate of the theory on its way to strong coupling is quite transpar- 
ent from the bosonized representation. First of all, the effect of the term 
STr M^^Mi with large couphng 92 is to "lock" the fields, since the function 
STr M^^M; has an absolute minimum at M^^M^ = 1; recall expression (pO]). 
Second, in the locked configuration M^ = Mi = M the terms multiplying 91 
and ^3 become r^(iV - 1) (^iSTr(MS3)2 + ^3STr2(ME3)). By parametriz- 
ing M = Te^T, with Y = +Ti^YJ]^ and S3TS3 = T~^, we can write them 
as 

^iSTr(MS3)2 + 53STr2(MS3) =^iSTre2^ + (?3 (STre^Sg)^ 

This is a potential for Y which, from (^), has an absolute minimum at 
Y = 0. (Here we again benefit from our careful construction of the target 
of the WZW model as a Riemannian symmetric superspace.) At strong 
coupling, we may set Y = 0. Doing so, and inserting Mi = T^ into the 
expression for the bosonized action 6*2, we obtain 

S'^ = NW[T^] + 9^ + 33 + N94 f d^^sTrdT^E^dT'E^ . 

(2vr)^ J 

We finally set q = TJ^^T'^ and identify q with the field of the nonlinear sigma 
model. Then, by a slight extension of the calculation of Section |^, we find 
that 5*2 reduces to the effective action Scs[q] given in ([2T|), with couplings 

f = {N-{9^ + 9s + N9,)/n)-' , 9 = ±Nn . (48) 

What will be the fate under renormalization of the nonlinear sigma model 
so obtained? Recall that for the choice of bare couplings we made, 91, 92 and 
93 flow to positive values whereas 94 becomes negative. By linearly combining 
the basic flow equations (|l^) we flnd 

^1 + ^3 + Ng, = JJ1±^ - (^i±^ _ (1 _ 2/N)9l - 2(9. + 92)93 ■ 
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The right-hand side of this equation is negative, so the combination gi + g^ + 
Ng^ has a negative rate of change, for all A^ > 2. (This happens in spite 
of the fact that the individual couplings increase in magnitude.) Hence, the 
nonlinear sigma model coupling in ( ^8]) is roughly of order 1/A^ and decreases 
under renormalization. It is therefore reasonable to expect that the model 
lies in the metallic phase, where the RG flow is attracted by the Gaussian 
fixed point / = 0. Our argument then says that the RG flow takes the 
Lagrangian L2 with bare couplings gi = (?2 and gs = gA = into the metallic 
phase. This remains true for arbitrarily small gi = g2, because this coupling 
is (marginally) relevant. Thus, no matter how close to the free-fermion point 
we start, the RG flow will go to the metallic fixed point. (Recall that this is 
not true for N = 2 with equal velocities, in which case the flow is attracted 
toward symmetry class All. However, we expect a finite anisotropy in the 
velocities to stabilize the flow inside class D.) In other words, for A^ > 2 we 
are sure of the existence of a metallic phase, and the free-fermion point sits 
right on the boundary of that phase, as claimed. 

On the other hand, the situation for A^ = 1 is qualitatively different. In 
that case, the only type of local disorder available in class D is randomness 
in the mass. If we set the statistical average of the random mass to zero, 
mo = {m{x)) = 0, there remains only a single coupling gM, which is defined 
by (HI). The bosonized action then reads simply 



S', = W[M] - -^ f (fx STr JS3 JS, 



This is a WZW model perturbed by a current-current interaction. Comput- 
ing the beta function directly from the OPEs for the WZW currents J and 
J, one finds that the perturbation is marginally irrelevant in the physical 
range gu > 0, in agreement with the flow equation (pISl). (The perturbation 
would be marginally relevant for negative g^j- However, for that sign of the 
coupling the C\D WZW model is unstable and does not exist; cf. the discus- 
sion at the very end of Subsection fr2| .) Thus, the RG flow is attracted by 
the WZW fixed point and the system in this sense is critical. 

The discussion so far assumed ttlq = 0. When a nonzero average mass is 
introduced, the bosonized action S'2 acquires an additional term 

{inio/i) /t/^xSTrMSg . 



This term is relevant by power counting at the WZW fixed point, and puts 
the system off criticality. Physically speaking, a finite Dirac mass rriQ opens 
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a gap around E = in the energy spectrum of the pure system, which causes 
locahzation of all low-energy states due to disorder (at least as long as qm is 
small enough). Therefore, the critical line segment mo = 0, qm > separates 
two insulating phases. As is known from [^, the distinction between the 



two phases is of a topological nature. If one of the two phases, say niQ < 0, is 
a plain insulator, the other one (mo > 0) can be likened to a quantum Hall 
fluid, in the sense that the presence of a boundary gives rise to chiral edge 
excitations leading to a quantized Hall-type response. (In the language of 
the 2d Ising model with weakly disordered bond strengths, the two phases 
correspond to the paramagnetic and the ferromagnetic phase.) Thus we have 
two insulating phases, and there is no room for a metallic phase, at least not 
in the vicinity of the free-fermion point, for A^ = 1. On grounds of continuity, 
the local density of states at zero energy in these phases is expected to vanish, 
and we believe the local Z2 degree of freedom of the C\D WZW model to 
play a crucial role in reproducing this feature in the field-theoretic formalism. 

By combining the information given, we arrive at Figure |l|, which draws 
a schematic picture of the phase diagram close to the free-fermion point. 
Three phases are seen to meet there: insulator, quantum Hall fluid, and 
metal. (Recall that in the physical application to superconductors, the ter- 
minology "metal" and "insulator" refers to i/ierma/ transport.) The two insu- 
lating phases are separated by a critical line, described by the field-theoretic 
model of Dirac fermions with random mass, or alternatively, upon nonabelian 
bosonization, by the C\D WZWfc=i model with a current- current interaction. 
The flow along the critical line terminates at the free-fermion point, which 
controls the critical behavior across the transition from the insulator to the 
quantum Hall fluid. The metallic phase exists for N > 2 species, or for N = 2 
with anisotropy in the velocities, but is absent for A^ = 1. The field-theoretic 
model for this phase is the CI | Dili nonlinear sigma model with couplings / 
and 9. It is expected |jTl| that the Hall response in this phase is not quantized 
but varies continuously with the angle 6. 

Two additional remarks are in order. First, a structurally similar phase 
diagram was suggested by Senthil and Fisher ||TT]|. The main difference is 
that these authors did not identify the multicritical point of class D as free 
fermions and, consequently, had part of the flow reversed. Second, recall 
from Section ^ that the topological term of the CI|Z^III nonlinear sigma 
model is trivial (at least in any continuum regularization of the field theory) 
for the case of one replica {n = 1), which contains all information about a 
single Green function at any energy E, and about a product of two Green 

58 



mass m o 



METAL 



CI I Dili 
nonlinear sigma model 





assless Dirac particle 



DiraC particle with random mass 
D WZW|j_, moderwith Jl Jr perturbation 




Figure 1: Schematic structure of the phase diagram close to the free-fermion point for class D 
in two dimensions. Field-theoretic realizations for some parts of the phase diagram are indicated. 
The coordinate along the horizontal axis is an A* -dependent function of the marginal couplings 
gi, ..., 5'4- For A* ^ 1, only the half-axis on the right of the free-fermion point is physically 
realized, with the coordinate being (^M ^ 9l ~^ 93 ~^ 94, ^0. The metallic phase exists for 
more than two species, or for N ^ 2 with anisotropy. 

functions at E = 0. It is hard to see how a term which is topologically trivial 
could drive a topological phase transition. This leaves two options: either 
the nonlinear sigma model (unlike Pruisken's model for the integer quantum 
Hall effect) is not a valid description of the critical line separating the two 
insulating phases, or else the driving agent for the transition is the Z2 degree 
of freedom of the model. In either case, the sole function of the topological 
term is to give rise to the edge current that is expected to exist for a system 
with boundary in the metallic and quantum Hall fluid phases. 



11 Vortices 

This is not yet the end of the story. It has recently been argued by Read and 
Green |]T3[ that vortex disorder has a significant effect on the phenomenology 
of quasiparticle transport and localization in class D. Without repeating the 
detailed discussion of |T^, we recall the following essential points. 

Basic to our study is the spinor field ip{x),il!^{x) of a Majorana fermion 
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with Hamiltonian 

H = (fx (^ipHdip^ + ipidip + mip^^^ 



which we here imagine to be the low-energy approximation to a mean-field 
Hamiltonian governing the time evolution of the quasiparticles in a chiral 
p-wave superconductor. Vortices are introduced by postulating a change of 
boundary condition: with (half-quantum) vortices present, it is no longer true 
that the Majorana spinor is a single- valued function of position. Rather, the 
spinor now reverses its sign on circling once around any one of the vortex 
singularities. (We assume the London limit where vortex cores are point-like 
objects.) The phase twist originates from gauging away the phase of the 
superconducting order parameter, which winds by 27r on going once around 
a vortex carrying half a magnetic fiux quantum. 

Thus the insertion of vortices amounts to the introduction of square- 
root singularities into the continuum Majorana field. Recall that the one- 
component Majorana theory in the massless limit (m = 0) is equivalent to the 
critical 2d Ising model. Under this equivalence, the square-root singularities 
in the Majorana field correspond to spin fields (or Kadanoff-Ceba disorder 
fields in a dual picture) placed in the partition sum of the Ising model. By 



using the known |HB| conformal dimensions of the spin field at criticality, and 
transcribing them to the supersymmetric setting. Read and Green were able 
to conclude that the perturbation of the massless free Majorana theory by 
randomly placed vortices is strongly relevant. 

The question then is: what does the relevant RG fiow go to? Read 



and Green |T^ suggested that the fiow leads into a metallic phase described 



by a nonlinear sigma model without the local Z2 degree of freedom which 
we have identified as being characteristic of class D. We have no direct 
proof of this conjecture at present, as we do not know how to incorporate 
square-root singularities into the supersymmetric extension of the nonabelian 
bosonization scheme. We can, however, give a partial verification as follows. 
Throughout our discussion, we are going to assume that the random 
mass m{x) is a smooth function of x. As a warm-up, let us start from 
the homogeneous limit m{x) = m-o > and insert a local inhomogeneity 
with the shape of a disk-like region V where the mass switches to negative 
values. Thus, m(x) is positive outside of V, negative inside, and zero along 
the boundary dT>. The zero-mass contour &D is a kind of "inner edge". As 



was nicely explained in |25], moving a Majorana spinor at energy E = Q 
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once around dT) reproduces the spinor with a phase shift of vr. Hence, Bohr- 
Sommerfeld quantization gives a bound state with low energy E = tt/L 
(actually, in the first-quantized formulation, a pair of bound states with 
energies ±ii^), where L is the circumference of the contour. The localization 
length of the bound state transverse to the contour is inversely proportional 
to the mass rriQ. 

At the critical point niQ = 0, a randomly but smoothly varying mass 
function produces a percolating network of zero-mass contours, and thus a 
macroscopic number of low-energy states. To capture the quantum physics 
of this random system, we turn to a variant of the Chalker-Coddington net- 



work 1^, originally designed to model the quantum percolation transition 
between plateaus of the integer quantum Hall effect. The original model 
involves random U(l) phases on the links, in addition to scattering at the 
nodes, of a square network. Since we are interested in the limit of zero en- 



ergy, where the states propagating along inner edges do not [^ carry any 
phases, randomness in the link phases is forbidden in the present context. 
Thus we are led to study a network model without any random phases, just 
randomness in the probability for scattering to the right {pr) or left (pl). 
At the symmetric point of the model, where Pl = Pr = 1/2 on average, the 
continuum limit is known [^ to be a Dirac fermion (or, equivalently, two 
copies of a Majorana fermion) with random mass and rriQ = 0. Thus we have 
come full circle and are back to our starting point. 

Consider now the effect of adding vortices to the system. As we recall, 
going once around a vortex in the continuum formulation twists the phase 
of the spinor wavefunction by minus one. For the case of an isolated zero- 
mass contour enclosing one vortex, the phase shift by tt has the important 
consequence of giving rise to a single zero- energy bound state or fermionic 
zero mode [^. In what follows, we shall be concerned with the effect of 
vortices on the critical system. This is captured in the network model by 
introducing "frustrated" plaquettes, by which we mean that the product of 
the link phase factors along any loop encircling such a plaquette equals minus 
one. Thus, randomly placed vortices in the continuum theory translate into 
randomly placed frustrated plaquettes in the network model. To create an 
isolated frustrated plaquette, we insert a semi-infinite string of minus signs on 
links, terminating at the plaquette. To create a large amount of frustration, 
we insert many strings. The maximal amount of disorder, corresponding 
to a high density of vortices, is realized by taking the link phases to be 
independent and identically distributed random variables drawn from the 
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uniform distribution on Z2 = {±1}, which is the same as the orthogonal 
group in 1 dimension, 0(1). In this hmit we arrive at the Chalker-Coddington 
network model with a single channel per link and Z2 invariant link disorder. 
The simplest choice is to take the scattering probabilities to be nonran- 
dom. Then, by the color-flavor transformation developed in |5^ for U(iV) 



(which readily extends [51| to O(A^) for all N including A^ = 1) or, alterna- 
tively, by Read's second-quantized setup [^, the Z2 phase invariant network 
model at its left-right symmetric point readily maps on the CI | Dili nonlin- 



ear sigma model with action (^I]) , and couplings / = 2/7r, 6 = it. However, 
there is one important difference from before, which is that the target space 
now has just one connected component. This can be understood by antici- 
pating from Section |1 1 . 1| that s G Z2 acts on the identity component by the 
trivial representation \s\ = 1, and on the other component by the faithful 
representation s = ±1. The latter is wiped out by averaging over the uniform 
distribution on Z2. (It reappears when the distribution is made nonuniform.) 
Since the free-fermion point formally corresponds (by saddle-point approxi- 
mation in the limit of small qm for A^ = 1; see Section |^) to the nonlinear 
sigma model with / = 1, and since the absence of a disconnected component 
of the target space does not change the perturbative RG beta function, the 
model with a weakened coupling / = 2/7r < 1 is expected to flow to the 
metallic fixed point / = 0. This result is compatible with the expectation 
of Read and Green that the addition of vortex disorder at the free-fermion 
point results in metallic behavior. By computing the functional integral of 
the nonlinear sigma model in zero-mode approximation, we find that the 
density of states for a finite network in the ergodic regime is 

p{E) = V + \b{E) . (49) 

Note that this answer differs from the ergodic limit of the density of states 
for quasiparticles in class D; cf. Section ^ and equation (^) below. 

Thus, vortices change the low-energy density of states in the metallic 
phase. The change is relatively minor there, being visible only on the micro- 
scopic scale of the mean level spacing i>~^. In the insulating phases, however, 
vortices are expected to have a much more pronounced effect. We have seen 
that, when vortices are absent, the local density of states at zero energy 
vanishes in that case. In contrast, the insertion of a finite density of vor- 
tices gives rise, by the above reasoning (for a slowly varying random-mass 
function), to an extensive number of quasiparticle states at very low energy. 
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In an insulating phase these do not communicate over large distances, and 
therefore they do not repel, leaving a finite density of states at E = 0. In the 
field-theoretic setting, we speculate that this comes about by the mechanism 
we identified above: vortices suppress the local Z2 degree of freedom of the 
nonlinear sigma model. 

We believe these differences in local observables and their field-theoretic 
manifestation to be sufficiently significant to warrant a refinement in vocab- 
ulary. What we propose to say is that Majorana fermions subject to vortex 
disorder belong to a generic (hybrid) class BD. This nomenclature is moti- 
vated by the following considerations. 

11.1 Random- matrix limits 

The fundamental point we are going to make, here, is that the Lie alge- 
braic structure underlying the Majorana fermion provides us with more 
than one kind of universal level statistics, and this suggests a refinement 
of the symmetry classification scheme. We begin by recalling from |^ that 
the basic random-matrix model for class D is the Gaussian ensemble over 
Dn = so{2N), the Lie algebra of the orthogonal group in even dimension. 
A Hermitian matrix drawn at random from iD^ has A^ pairs of eigenvalues 
±E, with the universal large- iV limit of the density of states being 

, , sin(27rz/£') , , 

Pz,W = . + 4^. (50) 

Note that this function is positive and smooth at ii^ = 0. 

Another Lie algebra, regarded as mathematically distinct from so(2A^) 
(as the root systems differ in structure) is that of the orthogonal group in 
odd dimension, denoted by B^ = so{2N + 1) in Cartan's notation. The 
distinction becomes expecially evident in the level statistics: the density of 
states of the Gaussian ensemble over iB^ has the large- A^ limit 

Pb{E) = u - ^^-^ + S{E) . (51) 

Here we see a singular term S{E), reflecting the fact that an antisymmetric 
matrix in odd dimension has one, and generically only one, eigenvalue at zero. 
(To understand this fact by way of example, recall that every proper rotation 
of Euclidian 3-space has one invariant vector, namely the axis of rotation.) 
Repulsion from that level is the cause of the quadratic law Pb{E) ~ E"^ near 
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E = 0. We mention in passing that the above law has also been found, in 
recent mathematical work [^], to be the universal limit of the density of 



zeros of suitably chosen families of algebraic L-functions. 

What we have reviewed and illustrated, then, is the fact that the large 
family of "orthogonal" Lie algebras splits into two classes, conventionally 
denoted by D and B. The latter class was mentioned in but not pursued 
there for want of a good physical example. Recently, it has been proposed 



5^ that the universal random-matrix limit of the density of bound states 
of a disordered vortex in a chiral p-wave superconductor, is given by class 
B. In earlier work [^, ^ |56[, the same limit for a conventional s-wave 



superconductor had been shown to be governed by class D. 

What is the relevance of all this to systems in the thermodynamic limit, 
when an extensive number of vortices is present? To answer that question, 
notice first of all that the ergodic limit (^) of the density of states for the 
Z2 phase invariant network model coincides with the arithmetic mean of the 
random- matrix answers for the classes B and D. This is already a first 
indication that the network model should not be assigned to the symmetry 
class D (at least not in a narrow sense), but is better placed in a hybrid class, 
which we propose to call BD. We are going to elaborate below. 

For the moment, we shall explain why the coincidence of (^91) with the 
arithmetic mean ^{po + Pb) is not an accident. Let us interpret |^ the 



network model as a quantum dynamical system evolving in time by discrete 
steps of size At. The discrete time-evolution operator U = exp —iHAt is 
a product of two unitary operators, one of which is diagonal (in the link 
basis) containing all the random Z2 factors associated with the links, while 
the other one encodes the random tunneling at the nodes. The eigenvalues 
of U are written e~*^^*, where E is called the quasi-energy. Let the network 
now be finite, with an even number 2N of links. Then U, being a unitary 
matrix with real entries, is an element of the orthogonal group 0(2A^). Note 
that for the model with Z2 phase invariant disorder, the probabilities for U 
to lie in the proper and improper parts of 0(2A^) are equal. 

Thus the time-evolution operator of the network model is a (sparse) ran- 
dom matrix on 0(2A^), and in order to understand the ergodic low-energy 
limit of the density of quasi-energies E, one needs to work out the random- 
matrix level statistics for that group. With some linear algebra, one can 
show that the level statistics for the two Haar ensembles on the proper and 
improper components of 0(2A^) coincide exactly with those of the classes D 
and B, respectively. (A noteworthy feature here is that, since the no ureal 
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eigenvalues of an orthogonal matrix come in complex conjugate pairs, every 
U in the improper component of 0(2A^), where Detf/ = —1, must have an 
eigenvalue —1, which then is accompanied by a corresponding eigenvalue +1, 
giving the delta function in (0).) The general answer for the universal limit 
of the density of states is some linear combination of po and ps as given by 



(50) and (^Tj). When the weights assigned to the proper and improper com- 
ponents of 0(2A^) are equal, the answer equals the arithmetic mean, which 
reproduces the result (^Q]), obtained from the nonlinear sigma model. 

Of course, the disordered quasiparticle systems we are interested in are 
typically very far from the random-matrix limit. One might therefore think 
that a nomenclature borrowed from that limit is too crude to capture the 
general situation. We insist, however, that the terminology we propose is 
fully supported from field theory, by the following argument. 

The key is to understand how the expressions (^,^) arise from the zero- 
dimensional limit of the field-theoretic representation by a nonlinear sigma 
model. As we recall, the target space of the field theory for class D has two 
connected components, and the result (^) comes from summing over both 
of them with equal Boltzmann weights. By inspection, one finds that the 
expression for class B results from a twisted target, where the Boltzmann 
weight of the component not containing the identity element carries a minus 
sign. It is useful to let the group Z2 act on the target space, so that twisting 
the Boltzmann weights is the same as acting with the nontrivial element of 
Z2. Now recall that the insertion of vortices amounts to creating frustrated 
plaquettes in the network model. Frustration of plaquettes, in turn, amounts 
to the presence of (strings of) links with minus signs. Any odd number of 
these will transfer the time-evolution operator U from the proper to the im- 
proper part of the orthogonal group, and vice versa. The crucial proposition 
we are driving at is that Z2 vortex disorder acts as Z2 sign disorder on the 
target space of the nonlinear sigma model. (Indeed, as we saw, going from 
D to S in the random-matrix limit is like switching between the proper 
and improper components of 0(2A^).) This proposition, albeit motivated by 
random-matrix considerations, does not depend on the random-matrix limit 
and we are therefore confident to predict its general validity. Thus, we expect 
that a square-root singularity (due to a half-quantum vortex) at a point p 
in the continuum Majorana field acts as the nontrivial element of Z2 at p on 
the target space of the CI|DIII nonlinear sigma model (whenever that model 
is a valid description); and, hence, averaging over randomly placed vortices 
suppresses the second target space component of that model. 

65 



This, ultimately, is the field-theoretic reason why vortices change the 
physics of disordered Majorana fermions, and why we believe that the sub- 
class D (without vortices) ought to be kept distinct from the more generic 
hybrid class BD. 

12 Open questions 

An extensive summary of our main results was already given in the intro- 
duction and will not be repeated here. Instead, we content ourselves with 
pointing out a few open questions we find interesting. 

i) It has been conjectured |Q that the spectral statistics of the bound 



states of an isolated disordered vortex in a chiral p-wave superconductor, 
coincides with the random-matrix statistics of class B in the universal limit. 
It would be desirable to try and prove that conjecture, by adapting the 
treatment of [|^ to derive the nonlinear sigma model for class B from a 
microscopic model. 

ii) To elucidate the effect of vortex disorder on a system close to the free- 
fermion point, it would be useful to know how to correctly transform the 
vortex (or square-root) singularities in the Majorana field to the bosonized 
representation. Since the insertion of a square-root singularity changes the 
Neveu-Schwartz fermions of the radial quantization scheme into Ramond 
fermions, the question is how to modify the C\D WZWfc=i model to account 
for a degenerate (Ramond) vacuum. 

iii) Although our work sheds much light on the 2d metallic phases of the 
classes D and BD, the insulating phases remain poorly understood. There, 
we expect domain walls in the Z2 degree of freedom of the CI | Dili nonlinear 
sigma model to be important in determining the low-energy physics. 

iv) According to ^^, the 2d random-bond Ising model (with fiuctuating 



signs, not bond strengths) can be cast in the form of a Chalker-Coddington 
network model with frustration (or vortices) on adjacent plaquettes. It would 
be interesting to understand better how this model fits into the expanded 
landscape of class BD versus D. In particular, the nature of the multicritical 
Nishimori point remains to be clarified. 
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